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14.  ABSTRACT 

This  report  results  from  a  contract  tasking  P.  N.  Lebedev  Physical  Institute  as  follows:  The  Project 
includes  five  Tasks: 

1.  The  modification  of  the  magnetosphere  by  the  controllable  flow  of  fast  electrons  coming  upwards  from 
the  artificially  modified  ionosphere.  Theoretical  investigation  and  the  conception  of  the  experimental 
program. 

2.  The  formation  of  super-narrow  ionosphere  structures  in  ionosphere  artificially  modified  by  powerful  radio 
waves  in  double-resonance  conditions.  The  super-narrow  structures  serve  for  the  strong  scattering  of 
UHF  radio  waves.  Theoretical  investigation  and  the  conception  of  the  experimental  program. 

3.  The  generation  of  powerful  gamma  emission,  powerful  bipolar  radio  pulses  and  optic  emission  from  the 
high-altitude  thunder  discharges.  Theoretical  investigation,  comparison  with  the  existing  observations  and 
plans  of  future  experiments. 

4.  High-energy  particles  in  thunder  phenomena.  Experimental  and  theoretical  investigation  of  intensive 
gamma  emission,  generation  of  HF  and  VHF  radio  waves,  electric  field,  electric  field  variations. 

5.  Transformation  of  an  extensive  atmospheric  shower  in  the  thunder  electric  field:  intensive  growth  of 
electromagnetic  component,  namely  of  the  high-energy  electrons  number,  of  gamma  quanta  and 
positrons  number.  Space  enlargement  of  the  shower.  Experimental  and  theoretical  investigation. 

Expected  results: 

1 .  The  theory  of  the  effects  arising  in  the  magnetosphere  modified  by  the  controllable  flux  of  fast  electrons 
coming  upwards  from  the  artificially  modified  ionosphere.  The  evaluation  of  the  possibilities  to  observe 
these  effects. 
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2.  The  evaluation  of  the  possibility  to  form  the  elongated  super-narrow  structures  in  the  ionosphere 
artificially  modified  by  powerful  radio  waves  in  double-resonance  conditions.  The  evaluation  of  the  UHF 
radio  wave  scattering  off  these  structures. 

3.  The  theory  of  the  high-altitude  thunder  discharges  generating  of  the  powerful  gamma  emission,  powerful 
bipolar  pulses  of  radio  emission  and  the  optic  emission. 

4.  The  theory  of  the  runaway  breakdown  for  the  case  of  very  high  values  of  electric  field  observed  in  the 
lightning  stepped  leader  which  led  to  the  effective  generation  of  high-energy  particles  and  intensive  fluxes 
of  gamma  emission. 

5.  The  results  of  detailed  simultaneous  measurements  of  the  set  of  parameters  associated  with  the 
lightning.  In  particular,  the  value  of  the  electric  field  and  it's  variations,  the  radio  emission  in  HF  and  VHF 
diapasons,  the  gamma  emission.  All  measurements  will  be  realized  from  several  registration  points. 

6.  The  results  of  simultaneous  measurements  of  the  set  of  parameters  associated  with  the  extensive 
atmosphere  shower  (EAS)  passing  through  the  thundercloud.  In  particular,  the  characteristics  of  high- 
energy  electrons,  muons  and  gamma  emission,  the  radio  emission  in  HF  and  VHF  diapasons  and  the 
value  of  the  electric  field  and  it's  variations.  All  measurements  will  be  realized  from  several  registration 
points. 

7.  The  theoretical  interpretation  of  the  measurements  associated  with  the  EAS  passing  through  the 
thundercloud. 

1.  The  theory  of  the  ultra  low-frequency  MHD  ionosphere  waves  generation  was  developed.  These  waves 
could  be  generated  by  the  action  on  ionosphere  by  modulated  microwaves.  The  exciting  currents  serving 
as  a  source  of  Alfven  and  magnetoacoustic  waves  were  estimated.  The  power  of  the  radiated  waves  and 
the  effective  coefficient  of  radiation  were  calculated  as  a  function  of  the  modulation  frequency. 

2.  The  theory  of  the  generation  of  the  super  narrow  striations  aligned  along  the  Earth  magnetic  field  was 
developed.  These  structures  can  arise  in  ionosphere  under  the  action  of  powerful  radio  wave  having 
frequency  close  to  both  upper  hybrid  frequency  and  to  a  multiple  gyromagnetic  frequency.  The  developed 
theory  was  compared  to  the  results  of  HAARP  experiments. 

3.  The  theory  of  the  high  altitude  atmosphere  discharge  arising  under  the  combined  action  of  the  runaway 
breakdown  in  a  thunderstorm  electric  field  and  the  passing  of  an  extensive  atmosphere  shower  was 
developed. 

4.  The  theory  of  runaway  breakdown  (RB)  for  the  case  of  a  strong  electric  field  was  developed.  The 
dependence  of  the  fast  electron  distribution  function  and  their  number  on  the  electric  field  value  was 
calculated.  The  numerical  simulation  of  electron  generation  was  done  taking  into  account  space  non¬ 
uniformity  of  the  electric  field. 

5.  The  investigation  of  the  correlation  between  cosmic  rays  and  the  bursts  of  radio  and  gamma  emissions 
was  carried  out  at  the  Tien-Shan  Mountain  Cosmic  Ray  Station.  Both  electric  field  and  its  variation  were 
registered  as  well.  The  data  analysis  discovered  cases  of  short-time  bursts  having  the  complex 
composition:  electrons  accelerated  up  to  about  10-30  MeV,  gamma  quanta  of  low  (tens  keV)  and  high 
(hundred  keV)  energies.  The  time  correlation  between  short  bursts  of  emission  with  the  lightning 
discharges  registered  by  the  electromagnetic  trigger  was  established.  The  space  correlation  of  these 
bursts  with  the  presence  of  electrically  charged  clouds  nearby  the  the  experimental  facility.  The  composite 
energy  spectrum  of  short  gamma  bursts  was  measured. 

6.  The  short  time  bursts  of  electron  flow  amplified  by  thunderstorm  electric  field  were  registered  directly.  They 
are  observed  exclusively  during  the  active  phase  of  a  thunderstorm,  absolutely  being  absent  during  fine 
weather,  rain  (even  heavy)  or  a  snowfall.  In  some  cases  the  correlation  of  accelerated  electrons  bursts 
with  the  moment  of  the  extensive  atmosphere  showers  (EAS)  passing  was  observed.  The  events  of  such 
kind  could  be  treated  as  the  direct  observation  of  the  lightning  discharge  generation  by  the  EAS  particles. 

7.  The  event  of  the  gamma  background  fall  correlated  in  time  with  the  thunderstorm  discharge.  This  event 
could  be  treated  as  the  first  direct  observation  of  the  earlier  predicted  form  of  the  discharge  (RB-EAS 
discharge).  This  type  of  discharge  occurs  in  thunderstorm  atmosphere  when  EAS  passes  through  the 
thunderstorm  electric  field. 
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Brief  description  of  the  work  plan 

The  Project  includes  five  Tasks: 

1 .  The  modification  of  the  magnetosphere  by  the  controllable  flow  of  fast  electrons  coming 
upwards  from  the  artificially  modified  ionosphere.  Theoretical  investigation  and  the  conception 
of  the  experimental  program. 

2.  The  formation  of  super-narrow  ionosphere  structures  in  ionosphere  artificially  modified  by 
powerful  radio  waves  in  double-resonance  conditions.  The  super-narrow  structures  serve  for 
the  strong  scattering  of  UHF  radio  waves.  Theoretical  investigation  and  the  conception  of  the 
experimental  program. 

3.  The  generation  of  powerful  gamma  emission,  powerful  bipolar 

radio  pulses  and  optic  emission  from  the  high-altitude  thunder  discharges.  Theoretical 
investigation,  comparison  with  the  existing  observations  and  plans  of  future  experiments. 

4.  Fligh-energy  particles  in  thunder  phenomena.  Experimental  and  theoretical  investigation  of 
intensive  gamma  emission,  generation  of  HF  and  VHF  radio  waves,  electric  field,  electric  field 
variations. 

5.  Transformation  of  an  extensive  atmospheric  shower  in  the  thunder  electric  field:  intensive 
growth  of  electromagnetic  component,  namely  of  the  high-energy  electrons  number,  of  gamma 
quanta  and  positrons  number.  Space  enlargement  of  the  shower.  Experimental  and  theoretical 
investigation. 

Expected  results: 

1 .  The  theory  of  the  effects  arising  in  the  magnetosphere  modified  by  the  controllable  flux  of 
fast  electrons  coming  upwards  from  the  artificially  modified  ionosphere.  The  evaluation  of  the 
possibilities  to  observe  these  effects. 

2.  The  evaluation  of  the  possibility  to  form  the  elongated  super-narrow  structures  in  the 
ionosphere  artificially  modified  by  powerful  radio  waves  in  double-resonance  conditions.  The 
evaluation  of  the  UHF  radio  wave  scattering  off  these  structures. 

3.  The  theory  of  the  high-altitude  thunder  discharges  generating  of  the  powerful  gamma 
emission,  powerful  bipolar  pulses  of  radio  emission  and  the  optic  emission. 

4.  The  theory  of  the  runaway  breakdown  for  the  case  of  very  high  values  of  electric  field 
observed  in  the  lightning  stepped  leader  which  led  to  the  effective  generation  of  high-energy 
particles  and  intensive  fluxes  of  gamma  emission. 

5.  The  results  of  detailed  simultaneous  measurements  of  the  set  of  parameters  associated  with 
the  lightning.  In  particular,  the  value  of  the  electric  field  and  it's  variations,  the  radio  emission 
in  HF  and  VHF  diapasons,  the  gamma  emission.  All  measurements  will  be  realized  from 
several  registration  points. 

6.  The  results  of  simultaneous  measurements  of  the  set  of  parameters  associated  with  the 
extensive  atmosphere  shower  (EAS)  passing  through  the  thundercloud.  In  particular,  the 
characteristics  of  high-energy  electrons,  muons  and  gamma  emission,  the  radio  emission  in  HF 
and  VHF  diapasons  and  the  value  of  the  electric  field  and  it's  variations.  All  measurements  will 
be  realized  from  several  registration  points. 

7.  The  theoretical  interpretation  of  the  measurements  associated  with  the  EAS  passing  through 
the  thundercloud. 

Technical  Approach  and  Methodology 

Tasks  1-3 

The  theoretical  approach  to  these  Tasks  is  based  on: 

1 .  Method  of  the  nonlinear  Boltzman  kinetic  equation. 

2.  Nonlinear  kinetic  theory  of  plasma. 
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3.  Analytical  analysis  of  the  nonlinear  equations  using  expansion  over  small 
parameters. 

4.  Numerical  simulations. 

Tasks  4  and  5 

1 .  Modification  of  the  Nal  scintillate  system.  Decrease  of  the  digitations  time  interval 
of  the  signal  registration  up  to  1  ps.  Extension  of  the  system. 

2.  Modification  of  the  registration  system  for  ’’slow”  and  ’’fast”  components  of  the 
electric  field  in  order  to  registry  them  with  10-100  ns  accuracy. 

3.  Extension  of  the  VHF  registration  systems  and  expansion  of  VHF  receiver 
bandwidth  up  to  1-3  MHz. 

5.  Expansion  of  the  Geiger-Muller  registration  system  up  to  effective  area  1  km2. 

6.  Creation  of  the  hadron  -  muon  detector  with  effective  area  up  to  0.5  km2  for  the 
fixation  of  EAS  coming. 

7.  Calibration  of  all  detectors,  in  particular,  using  standard  Csl37  sources. 

8.  Investigation  of  the  effect  of  the  thunder  electric  field  on  the  registration  system  in 
order  to  minimize  the  noise  and  to  discard  the  false  signals. 

9.  Registration  of  EAS,  gamma  and  radio  emission  in  thunder  conditions  and 
investigation  of  their  interconnection. 

10.  Investigation  of  the  interconnection  of  EAS  with  lightning  generation. 

Results 


Task  1:  Modification  of  the  magnetosphere  by  the  controllable  flow  of  fast  electrons 
coming  upwards  from  the  artificially  modified  ionosphere.  Theoretical  investigation  and 
the  conception  of  the  experimental  program 


The  small-scale  perturbations  of  ionosphere  plasma  arising  under  the  action  of  the  high- 
frequency  powerful  radio  wave  are  one  of  the  most  significant  phenomena  discovered  in  the 
modification  experiments.  These  perturbations  are  the  plasma  density  depletion.  A  large 
number  of  of  physical  effects,  from  anomalous  radio  wave  absorption  up  to  SEE,  are  connected 
with  these  perturbations.  An  ensemble  of  such  perturbations  creates  the  average  plasma  density 
depletion  aligned  along  the  magnetic  field.  The  investigation  of  the  effects  connected  with  such 
regions  of  plasma  density  decreasing  are  considered  in  the  Project.  Point  is  that  balance  of  the 
system  is  provided  by  the  currents  flowing  along  the  boundaries  of  the  regions.  If  the 
amplitude  of  the  pump  wave  is  modulated  with  some  frequency  the  modulation  of  the  plasma 
heating  and,  correspondingly,  the  oscillating  current  arise  in  turn.  These  currents  can  lead  to 
the  generation  of  MHD  waves  of  ultra  low  frequency. 

A  theoretical  investigation  of  the  excitation  of  MHD  plasma  waves  in  Earth 
magnetosphere  under  the  action  of  high-frequency  powerful  radio  waves  was  considered.  The 
high-frequency  wave  is  supposed  to  be  modulated  by  a  low  frequency  signal.  As  it  is  known 
the  powerful  radio  wave  form  striations  several  kilometers  long  in  the  ionosphere  plasma,  were 
the  electron  temperature  increases  for  2-4  times  [1].  The  heat  spreads  from  the  resonance  layer 
along  the  magnetic  field.  If  the  radio  wave  is  modulated  with  a  low  frequency  the  electric 
currents  are  excited  near  the  resonance  layer.  These  currents  could  serve  as  sources  of  MHD 
waves. 

We  assume  that  the  magnetic  field  is  directed  along  the  vertical  axis  z.  The  area  of  the 
radio  wave  acting  has  radius  R  ( R  is  about  10  km).  We  neglect  the  heat  transfer  across  the 
magnetic  field  and  the  recombination  losses.  The  resonance  layer  is  assumed  to  be  thin.  In 
these  assumptions  the  heat  transfer  equation  has  the  form: 
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I(ki£)  - b{T  -i]+ 5(z)(a+0ms{it) = %  (1) 

Here  r  =  T/Tx  ,  were  Too  -  the  temperature  out  of  the  heating  region.  The  longitudinal 
thermal-conductivity  coefficient  Kll  depends  on  the  temperature  as  «n  =  err5/2  [2],  the 
second  term  describes  the  energy  losses  in  collisions  of  electrons  with  neutral  molecules  and 
ions,  ,  b  =  8v  ,  were  ,  V  -  the  effective  collision  rate,  .  5  -  the  energy  fraction  lost  in 

the  collision  (it  is  about  me/mi ).  Third  term  describes  the  action  of  the  radio  wave,  Q- 
the  modulation  frequency  (it  is  about  1  Hz). 

Note,  that  due  to  plasma  quasi-neutrality  both  electrons  and  ions  take  part  in  the  diffusion 
(ambipolar  diffusion).  At  the  same  time  electrons  play  the  main  role  in  the  energy  transfer. 
That’s  way  the  temperature  perturbations  are  much  higher  that  the  density  ones,  which  we  have 
neglected  here. 

1)  Stationary  temperature  perturbation 

Let  us  treat  non-stationary  temperature  perturbation  as  a  small  one: 

t{z )  =  T0(2:)  +  A t(z.  t),  I  Ar|  <  T  (2) 

To  solve  the  equation  a  new  variable  is  defined 

V  =  z7/2,  y{z)  =  y0(z )  +  Ay(z,  t )  (3) 

It  leads  to  the  equation  for  Vo 

~Tj aV 0  +  KVo'  ~  1 )  =  a6(z)  (4) 

Using  natural  boundary  condition  Vo  ~ ^  at  •  Z/o  >  1  the  solution  of  (4)  should  be 
obtained 


ign(z)< 


a  (1 


b  \  9 


At  large  distances  the  perturbation  decrease  exponentially  with  the  characteristic  scale 

l  —  y/ajb 


The  condition  at  zero  should  be  also  taken  into  account: 


Lvw 

9y° 


2  _  7a2 
9  1 6ab 


,  Y0  =  yo{0) 


2)  Non-stationary  temperature  perturbation 


Equation  for  Ay  has  the  form 

„„5/7  d2Ay  ( 
dt  ay°  dz2  + 


=  -Y^7 (35(z)  cos  Qt 


Let  the  solution  form  to  be  [/  (.z) C 
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We  obtain 

r-4ff  =  ^o,  =  ^ 

vl'  a 

7  R 

/(+0)  _  /(_0)  =  -- 
Obviously,  the  asymptotic  behavior  of  the  solution  is: 


/(*) 


C 


oo,  Re  k  >  0 


(8) 

(9) 

(10) 


The  solution  of  the  equation  in  WKB  approximation: 

£ 

/  =  Cyl/We-"FM\  F(z)  =  j  - 


dr] 


vrw 


<ii) 


The  usability  condition  of  the  WKB  approximation  in  the  case  is:  |  at| />>  /  ,  in  another 
words  Q»b  (12) 

(Collision  frequency  in  the  ionosphere  F-layer  v  is  about  100  Hz,  8  is  about 

(30-1836M2-10'5,  thus .  5v  ~  10-3  ru .  <  Sly 


Constant,  C  is  determined  from  the  matching  condition  at  z=0: 

A/'(0)  =  C  ( 4y<rSS/28Ajfi(0)  -  2Ky0-5/28)  =  78 


28 

X 

c  = 


2  a 


7  0Yr 


5/28 


0 


|  aY0  9//14  +  4 atx 


=  \C\ 


Thus,  the  solution  for  k  has  the  form: 

A y{z)  t )  =  \C\y728e~KlF^z^  cos[Qt  +  ip  —  k2F(\z\ 


K  —  +  ili  2 


Finally, 


returning 


,  we  obtain  for  the  temperature: 

2  -5/7  A  2 


AtOM)  =  ^-0/'Ay  =  15/28e-~lF^cos[fi£  +  ^  -  «^F(|z| 


(13) 

(14) 

(15) 
to 

(16) 


The  characteristic  scale  of  this  distribution  is  ^  Vlftl  ^  \/a/-^  ~  km. 

In  reality  temperature  does  not  decrease  so  strong  out  of  the  heating  area.  There  is  some 
boundary  layer.  Its  width  could  be  estimated  using  the  fraction  of  longitudinal  and  transverse 

heating  transfer  «i/«i  =  kA ')  ,  were  -  the  electron  cyclotron  frequency.  For 

boundary  layer  width 
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we  have  "'ll/ 


K\\/L2  ~  K]_jh2  c  h  ~  L/(, 


/\^eTe).  in  our  case 


U)P  ™  107 


Hz,  thus  in  the  case 


Thus  we  have  estimated  the  temperature  oscillations  induced  by  the  absorption  of  the 
pump  mode  in  the  Earth  ionosphere.  Now  the  problem  of  the  oscillating  currents  should  be 
considered.  They  are  the  sources  of  magneto-hydrodynamic  waves. 

The  equation  of  continuity  and  the  equation  of  the  motion  in  the  frame  of  the  two  two- 
fluid  hydrodynamics  model  have  the  form  [2]: 

dne  _  , 

— — h  dw{neve)  —  0 

dt  v  e>  (17) 

drii 

— — h  awlriiVi)  =  0 

dt  v  l>  (18) 

dve  1  -*■  —*■ 

mene( +  (veS7)ve)  =  -Vpe  -  ene(E  +  ~[ve,  B ])  +  R 

at  c  (19) 

dv-  -  1 

miTii{ +  (viV)vi)  =  -Vp*  +  en^E  +  B})  -  R 

dt  c  (20) 

Here  we  have  neglected  the  viscosity  and  assume  Z  =  l. 

-  the  interaction  force  between  electrons  and  ions,  composed  of  the  friction  force  and  the 
thermal  force. 

Note,  the  Debye  radius  is  about  one  centimeter  in  ionosphere  plasma,  while  the  characteristic 
dimension  of  the  system  under  consideration  is  about  one  kilometer.  Thus  the  quasi-neutrality 
approximation  is  well  fulfilled  ne  =  ni  =  n . 

Let  us  turn  to  linear  combinations  of  the  obtained  equations  taking  into  account  the  relations 

H(  =  niandme«mi. 

Equation  for  hydrodynamic  velocity  f  =  {meve  +  'iThiVi)/(me  +  lfHi)  —  G 
(the  sum  of  equations  (19)  and  (20)): 

dV  „  1  r^jTi 

p—  —  — Vp  +  ~[jB\ 
dt  cl  1  (21) 


(here  F‘¥,  rFrFi, 

Generalized  Ohm  law  (the  difference  of  equations  (19)  and  (20),  previously  divided  on 
the  corresponding  masses): 

- £  - * 

E  +  t\VB]  +  —  (Vp,  -  flT)  =  -  +  —  +  — L IS] 

C  ^71  e  (7 1|  (J  j_  CTleC  (22) 

fj  I  =  Iff]  T  —  ]J  ^  (Jll  ”  2(7  I 

here  1  e/  e  (  6  ),  II  ;  thermal  force  for  the  case  of  strongly  magnetized 

plasma:  4  =  -0.71nV||7\ 

The  electric  field  is  potential  E  =  -Vtf  an(i  the  term  containing  gradient  of  the  pressure  leads  to 
the  renormalization  of  the  potential  only,  it  can  be  omitted. 

Note,  that  here  and  further  on  fields  E  ,  velocities  V ,  and  currents  J  represent  the 
fluctuating  parts  of  the  total  quantities  (i.e.  changing  proportionally  to  Just  the 
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fluctuating  parts  will  produce  MHD  waves  in  the  final  analysis).  The  term  should  be 
neglected  in  the  linear  approximation. 

Finally,  subtracting  (18)  from  (17)  and  taking  into  account  the  quasi-neutrality 
criterion,  the  current  closure  condition  could  be  obtained: 

div/=0.  (23) 

The  electric  field  is  determined  from  this  condition  in  the  case  of  quasi-neutral  plasma  [3]. 
Longitudinal  and  the  transverse  components  of  the  electric  field  (relative  to  magnetic  field)  are: 


->  ViiT^  ?'|| 

E\\+  0.7  11  =  J 11 

e  er,. 

(24) 

■  ~[VB\  =  ^  +  —  [j±B] 
c  <t_l  enc 

(25) 

using  V  ~  one  obtain 

_  iB2 

^  ^ 

(26) 

It  is  seen  from  the  equation  (21),  the  pressure  has  an  shock  at  the  heating  region  boundary.  It 

leads  to  the  origin  of  diamagnetic  current  7<p  circulating  along  the  cylinder  boundary  in  the 
layer  with  the  thickness  h  : 

.  rip  rip 

U  =  -^[VTeB]  =  ——T0At(z,  t)S(r  -  R)ev  ^ 

Turning  back  to  the  Ohm  law  and  substituting  (26)  into  (25)  we  obtain: 

E±_  +  Kh  =  j±/&±_  H — —\j±B] 

1  enc u  1  (28) 


V  -  — 


The  comparing 


Elm  . 


Ua  ~  ~  106 

mi 

~  -  ioo 

B  mi  1 2 


So,  the  main  role  belongs  to  the  resistance  .  Due  to  the  potentiality  and  axial  symmetry  of 

the  electric  field  ^  .  Thus  the  bulk  current  is  small:  3  <+>  ^  3r 

The  system  of  equations  is: 

^  0.7  dT  Jz 

Ez  H  v  — 

e  oz  a  (29) 

Er  77 jr  (30) 


The  equation  for  the  potential  should  be  obtained  using  ^  ~  and  the  closure  condition 
(23): 

oz  \  oz  J  77  r  V  or  J  e  c)z  \  oz  )  rm 
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Inside  the  heating  region  r  -  0.7T/e^ 

Outside  the  heating  region  we  have  the  equation: 

d  (  dv\  1  1  (  d(p' 


dz 


(•I) + a  (■ 


dr 


(32) 


where  (7  —  COflSt . 

The  coordinate  transform  T1  =  \fRffr  reduces  the  problem  to: 

A'<p= -4ttq(z),  (33) 

where  "  4re  ()£  ,  at  that  the  source  has  the  form  of  a  flat  “pancake”  with  the  thickness  Tv 
and  the  radius  . 

Let  us  make  use  of  electrostatic  analogy  —  the  problem  of  plane  capacitor  electric  field. 
We  are  interested  in  the  electric  field  behavior  at  its  border,  exactly  there  the  transverse 
component  does  arise,  and  thus  the  current  flows.  It  is  seen  that  if  we  do  the  inverse 


transformation,  that  currents  flow  in  the  thin  boundary  layer  with  thickness 


Li 


only. 


Let  us  consider  the  auxiliary  problem.  Let  the  charge  plane 
(  —  77<x<0,—  OO  <  y  <  +Oo)?  were  H  is  ver y  large,  have  the  plane  charge 

density  V  ? .  The  component  of  the  electric  field  Ex  in  the  point  with  coordinates  (^>^),  1>(I 
should  obtained. 

We  have: 

R 


a 


$  [-x  +  X)  dx  dy 


[-x  +  XY  +  i-y  +  Yy  +  i-z  +  Z?)^ 
/  It2 


(i  +  Ar)  dx 


X?  +  {Z-  zf 


=  ip  In 


x2  +  (z  -  zy 


For  the  neutral  is  average  layer  with  the  finite  thickness  and  charge  density  ,the  electric  field 
is: 


Ex  =  -  [  g(z)ln(X*  +  (z  -  Zf) 


dz 


(34) 


For  our  problem,  taking  into  account  that  during  the  back  transformation  electric  field  is 
transformed  as  dip) dr  =  y/TZadtp/dr  the  current  Jr  has  the  form: 


jr{r,  z)  =  J  In (7^.(7 Ar 2  +  (z1  -  zf)dz', 


where 


A  r  —  r  —  R. 


(35) 


Thus,  the  structure  of  the  currents  being  the  sources  of  electromagnetic  waves  was 
investigated.  Now  let  us  consider  the  problem  of  MHD  wave  emission. 

General  theory 

Let  us  consider  the  system  of  magnetic  hydrodynamic  equations 

dp  ^  ds 

¥  +  div(p,)  =  0,  -  =  0 


10 


(36) 
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dv  „  1  r  -r>  ^ 

Pdt=-Vp+^[l°lH’H] 


OH 

—  =  rot  [v,H\, 


div  H  =  0 


(37) 


(38) 


The  external  current  which  is  the  source  of  the  waves  should  be  singled  out  of  the  total 
current: 

rot  »  =  f(?+?)  (J9) 

Here  i  and  e  mark  the  internal  and  external  currents. 


Only  the  ponderomotive  force  should  be  included  in  the  equation.  It  affects  on  the  internal 
current,  what  corresponds  to  the  transformation  in  the  equation: 

^  47T 

rot H  — >  rot H - je 

c  (40) 

The  linearized  equations  could  be  obtained  by  considering  the  magnetic  field  and  the 
density  exiting  .  h  and  P  to  be  small  and  neglecting  the  quadratic  in  velocity  term 

,  {vV)v.  ~ 


dh  _  - 

—  =  rot[?7,  H],  div  h  =  0 


dp' 

dt 


+  p  div  v  —  0 


pdJt  =  ~U2°VP>  +  ^[r0t  K-'  ^ 


(41) 

(42) 

(43) 


In  the  case  the  sound  velocity  uo  y/  W/ dp  ^  a)  10  cm/s  js  jjmch  jess 

then  the  Al’fven  velocity-  di / \f^dp  ~  10  cm/s.  Thus  the  term  containing  the 

density  gradient  should  be  neglected  as  well. 


The  equation  on  .  V  could  be  obtained  by  differentiating  (43 


d2v  1  ,  1 

PW  =  4^  +  c 


H 


dj‘ 
’  dt. 


and  using  (41): 


(44) 


It  is  convenient  to  rewrite  (44)  in  cylindrical  coordinates,  where  it  decay  in  two  equations: 

d2vr  2f  &2vr  d2vr  1  dvr  1  \  H  djr 

- =  Vt  f  —  +  —  +  -  ^vr  ' 


dt 2 


hp 


(Jz2 

d2v 


dr 2  r  dr  rs 


pc  dt 


dt2 


=  V1 


d2v v  H  djr 


+ 


dz2  pc  dt 


(45) 


(46) 


The  wave  vector  ,  k  does  not  have  the  P  component  due  to  the  symmetry.  Thus  it  is 
seen  that  equation  (31)  describes  the  medium  oscillations  in  the  .  k  —  H  plane.  It 
corresponds  to  the  fast  magnetic  sonic  wave.  The  equation  (32)  describes  the  oscillations 
perpendicular  to  the  .  k  —  H  plane.  It  corresponds  to  the  Alfven’s  wave. 
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Alfven ’s  wave 

To  solve  the  equation  (46)  with  the  current  determined  by  (35)  the  Green  function  should  by 
found  as  the  solution  of  the  equation: 

d2g  J/2&g 
Adz 2 


-  V2AA  =  5{z)eint 


dt2  ~dzz  x  '  (47) 

Taking  into  account  that  the  required  wave  is  outgoing  one,  it  should  be  obtained  that: 

giflt— ik\z\ 

9  =  HikVf'  (48) 

where  ^  = 


The  solution  in  the  point 


Z  >  0,  Z  »  L 


vv  =  e 


iilt—ikZ 


ii 


2pcVA 


r  J  3r{z)t 


,ikz 


dz 


(49) 


Note,  that  •  VA/(l  ~  1000  km  I§>  L .  Thus  the  exponent  should  be  expanded. 

.2 


Taking  into  account  Q  ~  d~T j  dz  we  obtain 

r 

jr{z)  dz  —  jT(z)zdz  —  0, 


/ 


(50) 


r ) 

and  thus  the  first  nonzero  term  will  contain  Z" . 

The  wave  propagates  in  the  thin  boundary  layer  and  thus 


VS(r  FI).  Then  the 

right  part  of  the  equation  should  be  integrated  over  the  thickness.  As  a  result  omitting  the  factor 
exp  [iQl  ikz )  we  obtain: 


V  = 


0.7  H  k 


U2 


a 


8tt  pee  V\  V  77 


HI 


d~T 


dza 


- z 2  ln(7iaAr_  +  (zr  —  z)2)dz'  dr  dz 


(51) 


To  estimate  the  integral  we  suppose  —  Ar(0)e  I  d  ancj  use  qie  characteristic 


v'  -  L,  •  ‘~^r  ^  L I vlZcr  Then  the  velocity  should  be  estimated 

accurate  within  the  numerical  factor: 

3 

n  \  i~  /  \  /  fi  /  f  /  \  / 

V 


H  n2  ATuL3  A T0  /SUV  c 

~pAA/l  n  ~  eHL  \  vT  /  Va  J 


(52) 


Magnetic  sonic  wave 

Let  us  consider  the  equation  (45).  Let 


VT 


v .  The  nonuniformity  should  be  taken  into 


account  by  integrating  the  equation  in  the  vicinity  of  the  point .  r  —  R. 


3*v  =  vl 


dt 


&2vr  d2v 

+ 


dv 

dr 


dz 2 


\R+0 


1  dv 
fo2  +  rdr 


1 


:V 


dv 

dr 


r-  j 
A'a  til})  <9  At 


r  /  R 


(53) 


ft-0  — 


H2  dt  (54) 

The  Fourier  transform  over  .  0  and  ,  t  of  these  equations  should  be  done  so  that: 


V 


=  J  j  u{uj,k9,r)e* 


dkzdu 


(55) 
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(56) 


Equations  take  the  form: 


— u ru  —  Vi 


2  d2u  1  du 
—kzu  +  H — — - o 


9it, 

<9r 


l«+o 


dr 


ft-0 


dr 2 
47rnTo 


r  <9r  rx 

icuAr  =  U {to,  ks) 


(57) 


(58) 


.  2  ,2 
a  =  a; 


Marking  01  ^  f^A  and  substituting.  X  =  (XT  the  equation  we  reduce 

(57)  to  the  Bessel  equation  of  the  first  order: 

d2u  1  du  1 

+  7^  +  (l-T7)“  =  0 


or 


(59) 


77, 


do:2  x  dx 

The  solution  of  (57)  has  the  form 

E(oj.  kz)Ji(ar),  r  <  R 
R(uj,  kz)J\  (or)  +  jD(w,  (or),  r  >  R 

Substituting  this  solution  into  (58)  we  obtain: 

RJ^ocR)  =  BJi{aR)  +  DYi(aR) 

EJ[{aR)  +  U/a  =  BJ[(aR)  +  DY{{aR) 

These  are  two  equations  on  three  unknown  quantities  and  we  have  to  apply  an  additional 
restriction  on  coefficients  ^and  .  It  should  depend  on  the  sign  of  I1. 

1}  .  a2  >  0  ^  k2  >  UJ2/Vl 

The  solution  at,  r  >  R  must  give  the  outgoing  wave.  Quantity  can  be  taken  in  the  form 

,  u 

U  —  a{kz)S{uj  —  ft)  +  a*(kz)S(ui  —  ft)  (62) 

u  ~  H^Siu  -  ft)  +  h[2)S{cj  +  ft) 


(60) 


(61) 


Then  the  solution  should  be 


Coefficients  E,  D  are  also  taken  in  the  form 

D  D_S(u!  —  ft)  +  D+5(uJ  —  ft)5  and  then  the  additional  condition  could  be 
obtained: 

D-  =  iE-\  D+  =  —iE+  (63) 

Solving  the  obtained  set  of  equations  and  using  the  expression  for  Wronskian 
J\(x)Yi(x)  —  Ji(x)y\  (.t)  —c2/{t(X)  we  obtain  the  solution: 

,  j,  -ilRa{kz)Hf\aR)Ji{ar)5(ijj  -  Si)  +  k.  c.  •  +  Si),  r  <  R 

m  z,r)~  _  fi)  +  K.  c.  •  6{u>  +  Si),  r>R  (64) 

(k.c.  -  complex  conjugation) 


2)  •  a2  <0^k2>  uo2/Vl 

Let  P2  ~  -«2. 
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The  solution  finite  in  infinity  expressed  over  the  modified  Bessel  function  has  the  form: 

G(w,kz)h(fir),  r  <  R 
Q(lj,  kz)I<i  (fir),  r  >  R 


u(w ,  kz,  r ) 


(65) 


Substituting  it  into  (58)  we  obtain  the  set 

GIl(fiR)  =  QKl(fiR) 
GI((fiR)  +  U/fi  =  QK((fiR) 


(66) 


Finding  the  coefficients  and  using  the  expression  for  Wronskian 
Ii(x)K-±(x)  —  1 1 {x)K\  (x)  —  \/x^  we  obtain  finally: 

-Ra(kz)K\(fiR)Ii(8r)5(uj  -  Cl)  +  k.  c.  *  6(w  +  Cl),  r  <  R 

—Ra(kz)Ky (fir) I]  (fiR)S(uj  -  Q)  +  k.  c.  •  5(oj  +  CI ),  r  >  R 

(67) 


u(u),  kz,  r)  = 


We  found  the  asymptotic  expression  for  the  wave  at  large  ,3  ,  when  T  —  R  for 
simplicity.  According  to  (55),  (64),  (67): 

n/vA 

v(t,  z,  R)  =  2e 


.7T 
— l~r 


^  R  J  a(kz)Hi\aR)Ji(aR)cos(kzz)  dki 


Too 

-R  j'  a(kz)Ki(fiR)Ii(fiR)  cos(kzz)  dkz 


n/vA 


K.C. 


The  characteristic  value  of  the  argument  a(^z)  is  VF  Taking  into  account  that 
L  ~  R  Va/CI  ancj  fi  —  \kz\  we  can  supp0se  that  the  integrand  changes  just  a 
little  within  the  range  •  IT  Va/^]  we  obtain-  (x  ~  kzR,  7]  —  z/R ) 


V 


(t,  z ,  R)  =  2e  J a(x/R)Ki(x)Ii(x)  cos (xrfidx  —  k.c. 

o  (68) 


It  is  well  known  that  the  convolution  of  the  product  of  two  functions  is  equal  to  the  Fourier 
transform  of  these  functions.  The  Fourier  transform  of  a(K) 

was  calculated  earlier  (see 

(62)  and  (58)),  and  for  we  use  the  known  formula: 


Ki(x)ii(x)  cos(xy)dx  =  — - — — — - — 


(69) 


Obviously  y^oo 


The  function  At(?7)  concentrated  mainly  in  the  region  I V I  <  LfR  ~  1  When 
V  ^  1  we  have 
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~  iz(A-  A-  ‘t- _ 1\ 

(l'2'  Ar(V-y)dy~ 


|y|3 


-hoo 

/ 


+QO 


r  "  j  in- y) 

— OO  “OO 

It  is  seen  that  at  large  Z  •  v(t,r,z)  ~  \/z 


Resume 

1)  The  temperature  disturbances  of  the  ionosphere  arise  under  the  action  of  the  powerful 
radio  wave  modulated  with  low  frequency  ~  1-100  Hz.  They  propagate  as  waves  decaying 
on  the  length  about  L  ~  10  km. 

2)  Two  types  of  currents  appear  in  these  structures:  the  diamagnetic  current  circulating 
along  the  circle  near  the  tube  boundary  and  the  current  flowing  in  the  boundary  layer  in  the 
vertical  direction.  The  last  starts  at  the  middle  and  returns  on  the  tube  buts.  The  current  of  the 
first  type  excites  the  magnetic  sonic  wave,  while  of  the  second  type  -  the  Alfven  wave.  At  that 
the  Alfven  wave  is  channeled  inside  the  tube  while  the  magnetic  sonic  wave  propagates  in  all 
directions.  That  why  the  Alfen  wave  have  the  possibility  to  propagate  on  distances  larger  then 
the  magnetic  sonic  wave  in  spite  of  the  fact  that  its  amplitude  is  small  due  to  the  smallness  of 
transversal  conductivity  and  due  to  the  quadrupole  character  of  the  emission. 

3)  The  analysis  of  the  obtained  results  show  also  that  the  energy  fraction  transferred 
into  Alfen  wave  strongly  depends  on  the  modulation  frequency  increasing  with  its  growth.  On 
the  contrary,  the  amplitude  of  the  magneto-sonic  wave  practically  do  not  depend  on  the 
modulation  frequency. 
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Task  2.:  Formation  of  super-narrow  ionosphere  structures  in  ionosphere  artificially 
modified  by  powerful  radio  waves  in  double-resonance  conditions.  The  super-narrow 
structures  serve  for  the  strong  scattering  of  UHF  radio  waves.  Theoretical  investigation 
and  the  conception  of  the  experimental  program. 

Theory 

Theory  of  the  formation  of  super-narrow  structures  aligned  along  magnetic  field  was 
developed.  These  inhomogeneities  arise  in  ionosphere  under  the  action  of  powerful  pump  radio 
wave  when  its  frequency  is  close  both  to  the  upper-hybrid  frequency  and  to  the  multiple 
electron  gyro  frequency.  Upper-hybrid  waves  excited  in  these  conditions  are  trapped  in  the 
regions  of  local  plasma  density  depletions.  Moreover,  low-range  Bernstein  modes  are 
generated  as  a  result  of  the  decay  processes.  Bernstein  modes  being  phase-correlated  with 
stationary  upper-hybrid  waves  are  also  stationary.  Stationary  Bernstein  waves  produce  plasma 
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inhomogeneiteis  aligned  along  magnetic  field.  According  to  the  theory  these  ingomogeneities 
can  scatter  radio  waves  having  frequences  up  to  1  -  3  GHz. 

The  striations  appear  due  to  resonance  excitation  of  upper  hybrid  (UH)  waves  by  the 
powerful  pump  radio  wave  in  the  vicinity  of  upper  hybrid  resonance  point  Zuh  in  the  ionosphere 
determined  by  relation: 


co 


(i) 

2 

(CO  pe 


pe(^UH)  =  ~^e 

Here  coO  is  the  pump  wave  frequency,  cope  is  the  Langmuir  plasma  frequency 
Ane2Ne(z)/m)  and  Qe  is  electron  gyro  frequency  (Qe  =  eH/mc). 

Let  us  discuss  now  the  special  conditions  when  UH  resonance  point  Z?uh  is  close  to  the 
multiple  gyroresonance  nQe\ 

nQe,  o)2(ZdUH)  =  a2e(n2- 1),  n  =  3,4,5, 6, . 


CD, 


(2) 


Thus  in  the  vicinity  of  ^uh  double  resonance  conditions  are  fulfilled: 


co n 


co, 


In  the  vicinity  of  Zfuh  the  pattern  of  striation  is  changed  strongly.  First  of  all  according  to 
nonlinear  theory  [1]  in  resonance  region  the  number  density  of  striations  is  strongly 
diminishing  with  approaching  of  relative  pump  frequency  shift  Aco  =  co$  ~nQe  to  the  resonance 
point.  Exactly  at  resonance  Aco  =  0  striations  disappeared,  but  their  number  is  growing  rapidly 
at  small  shifts  Aco  =  (0. 005-0.  Ol^coo  to  the  both  sides  from  resonance.  Together  with  striations 
at  double  resonance  point  disappear  the  anomalous  absorption  both  pump  and  probe  waves 
propagating  in  the  resonance  region.  Naturally,  quite  analogous  behavior  demonstrates  the  field 
aligned  scattering  (FAS)  determined  by  striations.  The  experimental  observations  near  multiple 
gyroresonance  region  (MGR)  agree  well  with  the  theory  [2,3]. 

In  MGR  region  stimulated  electromagnetic  emission  (SEE)  is  strongly  changed  as  well. 
In  the  downshifted  region  Aco  <  0  only  downshifted  (DM)  structures  and  its  long  tail  in  SEE  is 
seen.  In  the  resonance  point  Aco  =  0  DM  structure  disappears  fully.  In  the  upshifted  region  both 
DM  structure  and  a  new  one  broad  upshifted  maximum  (BUM)  structures  are  seen.  The 
behavior  of  SEE  emission  in  MGR 

region  have  the  theoretical  explanation.  The  existence  of  BUM  feature  in  upshifted  region  is 
connected  with  the  excitation  of  Bernstein  waves. 

The  dispersion  equation  of  these  waves  according  to  the  theory  [4-5]  has  a  form: 
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Kn_x  (k2  p2)  =  0 


(3) 


Here  Qe  is  gyrofrequency,  p  is  electron  Larmour  radius,  k  is  perpendicular  wave  vector  and 
Kn(x)  ~  xn/(2nn\),  x«l  is  modified  Bessel  function. 

In  [6]  it  was  shown,  that  Bernstein  waves  with  k\\  »  k  can  be  excited  above  upper 
hybrid  resonance  layer  by  second  order  four- wave  interaction  process.  An  important  aspect  of 
the  BUM  is  that  the  frequency  of  excited  wave  is  higher  than  the  frequency  of  the  pump  wave. 
It  means  that  the  simple  three-wave  decay  process  cannot  be  responsible  for  this  effect.  That  is 
why  the  four  wave  interaction  was  supposed  already  in  the  first  quantitative  explanations  [7]. 
The  detailed  theory  was  developed  by  Huang  and  Kuo  [8].  This  process  involved  two  pump 
wave  “photons”,  upper  hybrid  plasmon,  an  electron  Bernstein  plasmon  and  driving  low 
frequency  oscillation  near  the  low  hybrid  wave  frequency.  The  second  harmonic  oscillation 
(2coq,  £o  ~  Q)  associated  with  the  HF  heater  wave  acts  as  the  pump  of  the  process.  The 
frequency  upshifted  upper  hybrid  wave  (002,  k^)  and  the  frequency  downshifted  electron 
Bernstein  wave  (co\,  k\)  are  the  sideband  and  the  decay  mode  of  the  decay  process, 
respectively.  A  low  frequency  electrostatic  oscillations  (cos,  k)  in  the  frequency  range  of  low 

y 

hybrid  waves  colh  =  (copeOe )  2  is  also  induced  by  nonlinear  coupling  with  HF  waves. 
Resonance  conditions 
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col  +  co)  =  co0  =  co2-  cos,  (2  cox  =  co  *  +  co2 ),  kl+k  =  0  =  k2-k  (4) 

are  supposed  to  be  fulfilled. 

This  process  is  considered  to  be  the  main  mechanism  for  generation  of  the  BUM  features 
observed  in  stimulated  electromagnetic  emission.  The  threshold  of  the  pump  electric  field  for 
this  parametric  process  has  a  form  [8]: 
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Here  coeb  is  electron  Bernstein  frequency. 

For  typical  ionospheric  conditions  in  the  vicinity  of  4th  gyroresonance  fo~  4fc  =  5-60  kHz 
and  for  f2  ~  /bum  =  20-80  kHz  the  threshold  field  Etb  =  0.2-0. 6  V/m,  the  corresponding 

instability  increment  y  ~  500-100  s-1.  It  means  that  elongated  perturbations  (k\\  <<k  ) 

associated  with  electron  Bernstein  could  be  really  excited  in  ionospheric  modification 
experiments. 

Let  us  emphasize  now  a  very  important  new  feature  of  the  effect.  Bernstein  waves  are 
excited  in  four-wave  interaction  process.  This  process  is  considered  in  the  theory  as  interaction 
of  free  waves  propagating  in  homogeneous  plasma  [2],  But  in  real  ionosphere  conditions 
situation  is  different.  The  upper  hybrid  waves  which  play  significant  role  in  four-wave  process 
are  not  free.  They  are  excited  due  to  resonance  instability,  trapped  inside  striations  and  cannot 
exist  outside  them  in  resonance  region  [1],  It  means  that  the  upper  hybrid  waves  exist  as 
standing  waves  inside  striations  only.  A  wave  moving  to  the  right  direction  has  a  strong  phase 
correlation  with  a  wave  moving  to  the  left. 

This  attribute  of  upper  hybrid  waves  affect  strongly  the  four-wave  process.  First — 
Bernstein  waves  could  be  excited  inside  striations  only,  where  UH  waves  exist.  Second,  the 
decay  conditions  (4)  for  upper  hybrid  wave  moving  to  the  right  and  to  the  left  direction  are 
identical.  It  means  that  the  Bernstein  waves  excited  in  decaying  process  of  two  pump  wave 
“photon”  into  upper  hybrid  plasmon  and  Bernstein  plasmon  are  standing  wave  as  well.  The 
excitation  of  standing  Bernstein  waves  is  very  important  because  this  standing  wave  can  serve 
as  a  source  of  new  quasi-stable  small  scale  density  perturbations.  Actually,  due  to  the  action  of 
striction  force  the  standing  waves  can  generate  small-scale  elongated  magnetic  field  density 
disturbances.  To  calculate  this  effect  let  us  consider  kinetic  equation  for  electrons  moving  in 
external  magnetic  field  and  in  standing  Bernstein  wave  E  =  Eeb  sin (kebx)  cos (coebt): 


df  df 

—  +  v  —  + 
dt  dr 
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Here  Qe  is  electron  gyrofrequency  and  coeb  is  frequency  of  Bernstein  wave.  The  last  term 
VN, 

*  i  Ar  is  the  ambipolar  input  of  ions,  where  N\  is  average  on  oscillations  plasma  density 

7Vo 

disturbances  and  7}  is  ion  temperature. 

To  calculate  the  striction  effects  it  is  necessary  to  solve  Eq.  (6)  taking  into  account 
terms  of  the  order  of  E2eb.  After  integrating  the  distribution  function  on  d3v ,  and  averaging  on 
oscillations  we  will  find  the  density  perturbations  M/M).  Fulfilling  this  program  one  can 
obtain: 


Nx  _  3 

N0~4m(Tl+Te) 
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Thus  we  see  that  field  aligned  density  perturbations  with  scale  of  Bernstein  wave  would 
be  excited  by  standing  wave  due  to  the  action  of  striction  force.  As  they  are  much  less  than  the 
striations,  it  is  naturally  to  call  them  super-small-scale  (SSS)  density  perturbations. 
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Let  us  estimate  now  the  amplitude  of  SSS  disturbances.  According  to  the  theory  the 
amplitude  of  Bernstein  waves  is  small  with  respect  to  upper  hybrid  waves  and  could  be 
estimated  near  /?-th  gyroresonance  as  [8]: 

Eeb  _  keb  {n^e  ~COeb) 


Jeb 

72 

^  uh 


kuh  (0Ji,h  ) 


Bearing  in  mind  that  due  to  resonance  pumping  the  amplitude  of  upper  hybrid  waves  is  much 
higher  than  the  pump  wave  amplitude  E0  and  taking  into  account  the  analysis  of  £„/,  fulfilled  in 

our  previous  work  [1],  one  can  estimate  the  amplitude  of  Bernstein  mode  Ea  ~  .  The 
estimations  of  threshold  field  amplitude  (5)  near  4-th  gyrofrequency,  for  example,  in  a 
frequency  shift  range  20-80  kHz  gives  Eq  ~  0.3-0. 4  V/m.  Thus  the  mean  value  of  SSS  density 
disturbances  (7)  could  be  roughly  estimated  as: 
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As  one  can  see  from  (7)  the  scale  of  these  disturbances  is  of  the  order  of  keb  ■  In  reality 
they  are  strongly  elongated  along  the  magnetic  field  and  have  cylindrical  symmetry.  In  Quarter 
II  we  have  considered  a  two  dimensional  model.  We  suppose  that  the  super-small  scale  (SSS) 
disturbances  are  located  randomly  inside  striations  and  have  profile 


8N(rL  )  <x  exp[-  (kebrL  )2  J 


The  drifts  of  perturbations  are  small  and  we  can  consider  SSS  as  stationary.  In  this  proposition 
the  fluctuation  spectrum  depending  on  perpendicular  to  the  Earth  magnetic  field  wave  vector 
k±  has  a  form: 
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Here  P  is  a  filling  factor  of  SSS  disturbances.  To  estimate  this  factor  it  is  necessary  to  remind 
that  SSS  disturbances  are  excited  and  packed  inside  striations,  the  scale  of  striation  across 
magnetic  field  ls  «  3-10  m.  Hence  the  number  of  small-scale  perturbations  is  of  the  order  of 
square  of  their  scales.  Thus  the  factor  P  could  be  estimated  as: 

The  value  of  perpendicular  wave  vector  keb  as  a  function  of  frequency  shift  V  =  ®o  ~  n^e 
from  n- th  resonance  could  be  obtained  using  dispersion  relation  (3): 


kebPe  = 


-i  i 


2  w!  ^ 
n2  -1  co 


pe , 


(2n-2) 


n  =  3,4,5,6, 


(10) 


To  estimate  the  value  of  keb  let  us  choose  the  typical  frequency  shift  =  0.0lcope  _  q^g 
fluctuation  spectrum  (9)  for  different  n  is  shown  in  the  Fig.2.1.  For  example,  the  value  for  n= 4 
is  kehpe  =  0.42. 
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Fig.2.1.  The  fluctuation  spectrum  ^ ^ 


as  a  function  of  k±pe  for  multiple  gyroresonances 


n=3,  4,  5,  6.  Below:  two  axes  in  GHz  for  different  electron  temperatures. 


The  field  aligned  scattering  at  SSS  perturbations  is  proportional  to  |$V/  .  One  can  indicate 

the  following  main  features  of  this  scattering. 

1.  Frequency  range 

Characteristic  scattering  frequency  for  n- th  multiple  gyrofrequency  is  fnc  =  c^et  ,  where  keb 
is  the  perpendicular  wave  vector  determined  according  to  (10).  Taking  a  typical  frequency  shift 
value  A  a,  =  0-01  a>pe  ?  we  find: 

fnc  =  3.3(A:x/?jGHz,  for  Te  =  1000°, 
fnc  =  2A{klPe )  GHz,  for  Te  =  2000° 

The  frequency  dependence  is  shown  at  the  Fig.2.1.  We  see  that  for  n=  4  -s-  6  resonances  the 
typical  scattering  range  is  up  to  2  4-  3  GHz. 


2.  Maximal  value  of  scattering  crossection. 

From  the  estimated  value  of  density  disturbances  (8)  one  can  find  the  density  fluctuations  as 
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Taking  into  account  that  for  usual  striations  'Tmx  ~  108  m2  [9,  10]  is  reached  at  (Sn  n)2  «  10  4 , 
we  estimate  the  SSS  field  aligned  maximal  crossection  as  °",nax  (SSS ) «  107  m2 


3.  Frequency  spectrum 

As  is  well  known  the  field  -  aligned  frequency  spectrum  width  for  a  wave  scattered  by  usual 
striation  is  extremely  narrow  ~  0.5  Hz.  That  demonstrates  a  high  stability  of  the  whole  striation 
structure  developed  in  ionospheric  modification. 

The  width  of  frequency  spectrum  for  the  SSS  scattering  could  be  different.  According  to 
nonlinear  theory  there  is  a  strong  heating  of  electrons  inside  striations  [1],  In  average: 


A Te  ~  0.57],.  n  means  that  one  can  expect  a  significant  drift  motion  of  SSS  fluctuations  in  the 
direction  perpendicular  to  magnetic  field  B: 
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Taking  the  typical  scales  of  striations  l  ~10  m  [1 1],  one  can  obtain  from  (12) 
cm/s  and  the  typical  scattering  frequency  width  is 


5  x  103 
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fo 
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Thus  the  frequency  spectrum  width  for  field  aligned  scattering  of  UHF  radio  waves  at  SSS 
few  hundred  Hz. 


(13) 
is  a 


4.  The  growth  and  decay  of  SSS  scattering 

Let  us  consider  the  scattering  process  in  the  case  of  sharp  put  on  or  put  off  powerful  pump 
wave. 

After  the  put  on  of  the  pump  wave  the  striations  and  SSS  are  growing  together.  It  means  that 
the  characteristic  growth  time  of  SSS  scattering  is  of  the  order  of  a  few  seconds  -  the  same  as 
the  scale  time  of  striation  growth. 

Quite  a  different  picture  is  for  decay  time.  Striations  disappear  as  usually  in  a  few  seconds,  but 
SSS  density  inhomogeneities  are  created  by  striction  forces  and  have  a  very  small  scales 

r±  ~  (3-5)  cm.  The  plasma  diffusion  in  such  scales  goes  fast:  the  fluctuations  of  electrons 
disappear  due  to  diffusion  along  inhomogeneities  and  ions  -  due  to  free  motion  across.  It  gives 
the  life-time 

l  r 

r,,  ~  -A—  «  5xl0“2  s,  r± ~  —  ~10“4s. 

I|e  2  De  Vt 

So  the  life-time  of  SSS  after  the  pump  put  off  is  only  50  ms  or  less. 

Comparison  with  the  experiment 

According  to  the  earlier  elaborated  theory  powerful  radio  waves  which  frequency  is  close  to 
the  multiple  gyro  resonance  being  used  to  effect  ionosphere  produce  super  narrow  striations 
elongated  along  the  Earth’s  magnetic  field  direction.  According  to  our  suggestion  the 
experiment  devoted  to  the  search  for  such  inhomogeneities  was  realized  at  HAARP  during 
2008  spring  campaign.  The  ionosphere  heating  was  done  under  the  conditions  corresponding  to 
the  third  gyro  resonance.  The  GPS  signal  from  the  satellite  transited  over  the  HAARP  in  the 
experiment  moment  was  used  as  a  diagnostic  one.  The  phase  shift  of  the  high-frequency  radio 
signal  was  observed.  The  data  are  processed  now. 

We  have  proceed  the  analysis  of  the  phase  shift  of  the  GPS  signal  from  the  satellite 
transiting  above  the  zone  exited  in  ionosphere  by  the  powerful  radio  wave  having  frequency 
close  to  the  third  gyro  resonance.  The  data  analysis  shows  that  the  phase  shift  sharply  changes 
and  even  oscillates  with  the  period  about  20  s  just  after  heating  switch  on.  This  phenomenon 
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was  observed  in  four  of  the  six  sets  of  observations.  The  electron  density  perturbation  has  a 
huge  value  -  about  20-40  %.  The  presence  of  such  strong  long  scale  (1-5  m)  perturbations  is 
unlikely  from  the  theoretical  point  of  view.  Apparently,  effect  is  connected  with  the  dispersion 
of  the  GPS  signal  on  the  super  small  scale  irregularities. 

The  processing  of  HAARP  experimental  data  was  finished.  The  experiment  was  devoted  to 
find  the  scattering  of  the  GPS  signal  by  the  electron  density  perturbations  in  the  F2  region  of 
the  ionosphere.  The  effect  was  earlier  predicted  in  the  Project.  The  experiments  were 
conducted  using  the  HAARP  heater  having  the  radiating  frequency  / which  matches  3/b,  i.e., 
triple  the  local  electron  gyro  frequency.  Such  frequency  is  expected  to  generate  super  small 
irregularities  of  the  electron  density  which  can  scatter  GPS  signals.  It  was  found  that  the 
differential  phase  of  the  probe  GPS  signals  changed  abruptly  in  about  10  s  after  the  start  of  the 
HF -heating,  and  then  oscillated  with  the  heating  period  20  s.  The  oscillations  lasted  for  4-5 
minutes  and  then  disappeared,  presumably  when  the  resonance  condition  f  =  3/b  was  not 
satisfied,  although  the  HF-heating  continued. 

The  smallest  scale  resolvable  in  the  GPS  phase  data  are  of  an  order  of  tens  of  meters  given  the 
scan  velocity  of  the  GPS  beam  in  the  F2  region  and  the  2  Hz  sampling  rate.  Thus  the  experiment 
does  not  provide  a  direct  proof  that  the  observed  GPS  perturbations  are  caused  by  super  small 
scale  (SSS)  irregularities  such  as  predicted  by  the  theory.  However,  there  are  some  indications 
that  SSS  irregularities  play  a  role  in  this  effect  as  discussed  below. 

In  fact,  let  us  estimate  the  phase  change  due  to  GPS  signals  of  frequency  co  passing  through  the 
perturbed  region  of  the  ionosphere: 

.  ,  CO  f  ,  co  co]  8n 

A  <p  =  —  \£dzts - — — 1 

c  J  c  co  n 

2  V2 

where  s  is  the  refraction  index,  coe  =  (4  7ie  n/m)  is  the  electron  plasma  frequency,  Sn  is  the 
electron  density  perturbations  due  to  the  HF-heating,  and  I  is  the  l^nath  of  the  nprtm-hed  region. 

For  the  heating  frequency  equal  to  4.3  MHz,  and  assuming  that  coe  »2ttx  4.3  x  10  s  it  could 
be  obtained  that 

Sn  ~  10A^(raJ) 
n  I(km) 

The  differential  phase  observed  in  the  experiment  was  about  0. 1  radians,  while  the  width  of  the 
disturbed  region  does  not  exceed  5-10  km.  This  gives  the  electron  density  perturbations  of  10  - 
20%  which  is  almost  an  order  of  magnitude  higher  than  that  usually  observed.  Such  enormous 
phase  perturbations  occurring  on  a  one-second  timescale  indicate  that  scattering  of  GPS  signals 
by  the  SSS  irregularities  could  contribute  to  the  differential  phase  observed. 

Beside  the  abrupt  relaxation  of  STEC  perturbations  detected  in  experiments  (see  the  pulses 
marked  by  the  arrows  on  Figure  2.2  (bottom)),  shows  that  STEC  relaxes  to  its  unperturbed 
value  within  a  fraction  of  a  second  after  the  termination  of  the  HF  pulse.  The  fast  relaxation 
implies  that  it  is  caused  by  the  ponderomotive  force  associated  with  the  SSS  irregularities,  since 
the  relaxation  of  the  large  scale  irregularities  occur  on  a  much  longer  timescale.  Apparently  no 
other  process  in  the  ionosphere  can  lead  to  such  fast  scattering  of  GPS  signals. 

Thus,  the  observed  phase  oscillations  indirectly  indicate  the  presence  of  super  small  scale 
irregularities  of  electron  density  produced  by  the  HF  heating  at  the  frequency  corresponding  to 
3/b- 
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Figure  2.2  (top)  Absolute  STEC  versus  the  observational  time  is  revealed.  The  encircled  area 
shows  oscillations  with  20  s  period,  (middle)  The  5  minutes  period  of  STEC  oscillations  is 
zoomed,  while  (bottom)  the  de-trended  STEC  during  the  same  period  is  shown.  The  beginning 
and  end  of  the  HF-heating  are  marked  by  the  vertical  dashed  lines.  The  black  horizontal  bars 

show  the  heating  periods,  each  is  10  s  long. 
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Task  3:  Generation  of  powerful  gamma  emission,  powerful  bipolar  radio  pulses  and  optic 
emission  from  the  high-altitude  thunder  discharges.  Theoretical  investigation, 
comparison  with  the  existing  observations  and  plans  of  future  experiments 

The  theory  of  the  high-altitude  atmosphere  discharge  was  developed  conformably  to  the 
combined  action  of  the  runaway  breakdown  (RB)  and  extensive  atmosphere  shower  (EAS)  in 
thunderstorm  electric  field.  EAS  can  be  initiated  at  the  heights  13-20  km  only  by  particles 
propagating  with  small  angle  to  the  horizon.  EAS  weight  can  reach  one  kilometer.  It  was 
shown  in  the  developed  model  that  runaway  breakdown  is  essentially  amplified  under  these 
conditions. 

Introduction 

Runaway  breakdown  in  air  (RB)  is  stimulated  by  the  presence  of  high  energy  cosmic 
ray  secondaries.  Extensive  atmospheric  showers  (EAS)  are  accompanied  by  an  effective  local 
growth  of  the  number  of  cosmic  ray  secondaries  and  thus  have  a  strong  influence  on  the  RB 
process  [1],  The  combined  action  of  runaway  breakdown  and  EAS  lead  to  the  development  of 
RB-EAS  discharge  -  the  discharge  where  relativistic  electrons  play  a  decisive  role.  RB-EAS 
discharge  is  accompanied  by  strong  exponential  growth  of  the  number  of  energetic  and  thermal 
electrons,  positrons  and  gamma  quanta  [2],  It  can  serve  for  the  generation  of  a  strong  bipolar 
radio  pulse  [3]. 

A  distinct  class  of  radio  pulses  generated  in  thunderstorms  was  effectively  studied 
during  recent  years  in  observations  of  Smith  et  al  [4-6].  These  studies  allowed  establishing  that 
the  radio  pulses  having  enormous  peak  power  up  to  100  GW  are  emitted  in  the  wide  frequency 
range  by  the  intracloud  discharges  in  the  tropopause.  The  pulses  are  short  time  (<  10  ps)  and 
have  a  definite  bipolar  form.  That  is  why  they  were  called  narrow  bipolar  pulses  (NBP).  The 
observations  show  that  these  strong  radio  pulses  are  isolated  -  not  accompanied  by  usual 
lightning  leader  and  return  stroke.  Their  optic  emission  is  very  weak.  A  detailed  analysis  of  the 
whole  complex  of  observational  data  allowed  Jacobson  (2003)  to  state  that  NBP  is  a  new  type 
of  thunderstorm  discharge  quite  different  from  usual  lightning  [7],  He  speculated  that  it  could 
have  relevance  to  runaway  breakdown  effect. 

A  theory  of  RB-EAS  discharge  in  air  was  compared  with  the  results  of  NBP 
observations  and  qualitative  agreement  was  demonstrated  [8].  That  allowed  to  state  that  RB- 
EAS  discharge  can  serve  as  a  background  of  NBP  phenomena.  Note,  that  in  [8]  a  simplified 
model  where  a  high  energy  cosmic  ray  particle  moving  along  the  vertical  direction  parallel  to 
the  thundercloud  electric  field  E  was  used. 

In  the  same  time  the  NBP  are  generated  at  tropopause  heights  (10-20)  km.  At  such  a 
high  heights  energetic  cosmic  ray  (CR)  particle  has  a  well  developed  extensive  atmospheric 
shower  (EAS)  only  if  its  momentum  has  a  large  angle  6  to  the  vertical  direction  0  >  60°  -  70° 
[9],  The  first  goal  of  our  investigation  is  to  develop  a  theory  which  takes  this  fact  into  account. 
We  will  demonstrate  that  inclination  of  CR  particle  trajectory  lead  to  elongation  of  the  path  of 
runaway  breakdown  what  serves  for  significant  amplification  of  RB.  It  means  that  the  RB-EAS 
discharge  current  grows  effectively  with  the  inclination  angle  6  and  in  NBP  conditions  can 
reach  extremely  high  values  -  up  to  100  kA  and  even  more. 

Note  that  in  the  considered  previously  RB-EAS  discharge  models  the  electric  field  in 
thundercloud  was  supposed  to  be  a  given  stationary  function  of  height  z  [2,  8],  That  is  correct 
for  not  too  strong  current  only.  Let  us  remind  that  the  energy  spent  during  RB  -  EAS  discharge 
comes  from  the  electrostatic  energy  of  thundercloud  electric  field.  This  is  the  only  energy 
source.  It  is  evident  that  if  the  discharge  is  intensive  enough  a  part  of  electrostatic  energy  is 
dissipated  and  the  electric  field  in  the  cloud  should  diminish.  This  will  affect  the  RB  -EAS 
discharge  and  determine  its  nonlinear  modification. 
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The  second  goal  of  the  investigation  is  to  develop  the  RB-EAS  discharge  theory  taking 
into  account  nonlinear  process.  We  will  demonstrate  that  high  power  discharges  are  strongly 
modified  due  to  this  nonlinearity. 


RB  -  EAS  discharge  at  tropopause  heights 

Extensive  atmospheric  shower  (EAS)  becomes  well  developed  when  the  energetic  cosmic  ray 
(CR)  particle  pass  a  significant  length  r  in  the  atmosphere  r>  r£\  12].  Here: 


I 

^,\p(2)dz 


is  the  integral  of  atmospheric  density  p(z)  measured  in  g/cm2  and  6  is  the  inclination  angle  of 
CR  particle  having  energy  s.  The  radiation  length  r£  is  growing  with  energy  logarithmically 


r  =  In 


80MeV 


and  for  s=  1016  -  1019  eV  takes  the  value  r£  =  18.5  -  25.4. 

Atmosphere  at  the  heights  of  NBP  generation  (10-20)  km  is  not  dense.  For  example  in 
the  case  of  vertically  moving  CR  particle  for  the  height  z  =  16  km  the  integral  length  is  only  r 
=  5  g/cm2.  That  is  why  EAS  become  well  developed  only  if  the  energetic  cosmic  ray  particle 
momentum  is  directed  close  to  the  horizon  9h  =  90°  -  0  <  9*(e).  The  critical  angle  9*  for  the 
height  z  =  1 6  km  is  presented  in  the  Table  1 . 

_ _ _  Table  1.  EAS  critical  angles 


e(eV) 

9* degree 

0*rad 

1016 

16° 

0.28 

1017 

14° 

0.24 

1018 

12.5° 

0.22 

1019 

11.4° 

0.20 

Taking  the  electric  field  in  the  thundercloud  directed  as  usual  close  to  the  vertical  one  can  see 
that  now  the  RB  -  EAS  discharge  is  developing  in  quite  a  new  condition.  In  models  considered 
previously  the  CR  particle  momentum  was  always  supposed  to  be  directed  along  the  vertical 
axis  z  [2],  [8]  what  automatically  lead  to  the  development  of  the  discharge  along  electric  field 
Ez  with  the  velocity  of  light  c.  Now  the  CR  particle  is  moving  practically  perpendicular  to  E. 
How  will  the  RB-EAS  discharge  develop  in  these  conditions? 

First  we  have  to  remind  that  to  realize  RB  two  conditions  are  needed: 


-  a  seed  high  energy  (s>  0.5  Mev  )  electrons, 

-  the  electric  field  E  >  Ec  in  the  thundercloud. 

The  first  condition  is  fulfilled  for  EAS  because  it  generates  a  large  number  of  high  energy 
electrons  having  a  widely  spread  angular  distribution  in  energy  region  £  <  IMeV  [4,  12,  13]. 
Namely  this  energy  region  is  the  most  interesting  for  excitation  of  RB.  So  we  see  that  if  E  >  Ec 
the  both  requirements  are  fulfilled  and  RB-EAS  discharge  exists  in  the  considered  new 
conditions. 


But  there  is  a  specific  discharge  peculiar  property,  different  from  the  previous  case  v\\E. 
It  is  plotted  qualitatively  at  the  Fig.3.1.  In  Fig. 3.1  the  CR  particle  is  shown  moving  along 
x-axis.  Electric  field  is  directed  up  along  z  perpendicular  to  x.  The  RB  process  begins  from 
point  x  =  0  and  z  =  zc  where  E(x,z)  >  Ec  and  goes  down  along  z-axis.  It  is  shown  by  the  arrows. 
The  velocity  of  RB  region  boundary  in  z  direction  is  determined  by  distribution  function  of 
runaway  electrons  in  the  energy  region  £n0.3-0.5  MeV.  According  to  numerical  calculations 
Symbalisty  et  al  [12],  Lehtinen  et  al  [13]  the  average  directed  velocity  <vz>  in  this  region  is 

vz  =  2 7rjvpf(v,p,t)v2dvdp  ~  (0.3-0.5)c 


We  see  that  <vz>  is  significantly  less  than  c.  Due  to  this  as  the  RB  growth  time  rrb  depends  on 
E/Ec  only  (see  [14])  the  RB  effective  avalanche  length 
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laeff  T,-b  (2) 

is  diminishing.  This  diminishing  of  effective  avalanche  length  laeff  lead  to  exponential 
amplification  of  RB-EAS  discharge  in  a  given  thundercloud  electric  field  layer  L  (oc 
exp(L/laeff)).  Another  significant  amplification  process  is  connected  with  the  growth  of  EAS 
cross  section  which  is  determined  by  the  diminishing  of  electron  energy  losses  at  tropopause 
heights  due  to  low  density  of  air  molecules. 


Fig.3.1 
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The  kinetic  theory  of  gamma  emission  propagation  in  non-  uni  form  atmosphere 
was  developed.  The  goal  of  the  investigation  is  to  clear ify  the  nature  of  TGFs., 
gamma-bursts  registered  on  l  he  satellites  COMPTON  and  RIIESSI  as  comming 
from  the  Earth.  The  emission  was  found  to  be  assotiated  with  lightning,  A 
theoretical  model  of  gamma  radiation  propagation  and  its  yeild  from  the  non- 
uniform  Ear  til  atmosphere  is  needed  to  estimate  radiation  parameters  in  the 
source.  The  kinetic  equation  discribing  the  photon  propagation  was  obtained 
for  the  case  of  arbitrary  non  mini  form  media.  In  the  com  monease  it  has  the 
form: 

#n(r,  k.  t)  kdn(r.kA)  f  $n \ 

di  +Ck  dr 

where  n(r.  k,  i)  photon  density  in  the  phase  space  (r.  k).  (77)  ,ifl  f.  photon 
electron  collision  integral. 


(/(p')n(k')[l  +  ri(k)]  -  /(p)n{k)[i  +n(k')])  x 


xwfp'.k',  p.  k )tPptPpttPk\ 


(1) 


here  w(p'.  kf.  p.k)  the  probability  of  the  photon  transition  from  the  state  with 
impulses  (p.k)  into  the  state  with  impulses  (p',k').  In  the  relativistic  invariant 
form  this  probability  could  l>e  written  as: 


A  A 


e4  JUp  +  k-p1  -k1) 

^(p',  k'.  p.  k,  =  -R - — - 

_  /  m2  tn2\2  A  D 

fi  =  (1  +  X'"oJ  +B  +  a_1- 

£  electron  energy.  A  X  B  kinematic  invariants: 


A  =  kp  =  k(t  —  pcQ8%)  =  k'p'* 

B  =  k'p  =  kf(e  —  pco$£)  =  kp'. 

After  integrating  (1)  over  scattered  electron  impulses  pf ,  and  the  over  the  target 
poton  impulses  kf  the  collision  integral  takes  the  form: 


J  t{e  -  pcosQ 

where 


R  = 


^  m2(l—co$x) 

(*  —  peo$£)(t  —  pco$£f) 


+ 
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m2(l  —  co$£f)  e  —  pco$£'  +  k(l  —  co$x) 

e  —  pcoa^  +  k(l  —  coax)  m2(l  —  co$£f) 

X  -  the  angle  between  vectors  k  X  k',  sf  -  the  angle  between  vectors  p  X  k.  - 
tile  angle  between  vectors  p  X  kf,  R  —  R(k).  k  has  the  form 

r  _  k(e  -  pco&£) 

e  —  pcoa$f  —  k{  1  —  coax) 

In  the  case  when  lire  photon  energy  is  much  higher  then  the  electron  temperature 
the  electron  distribution  function  has  the  form: 

/( p.z)  =  N(z)6(p).  (3) 

Taking  (!$)  into  account  the  ex  pres  ion  (2)  could  be  rewritten  as 


2 m2k 


2^2  ^ (") J 


dfi'  coa2\ - (1  —  + 


1 


rn 


1- 


x 


x  ■ 


k2n(k.  n!) 


(1-  ^{1- coax))' 


9  k 

coa2x  H — -(1  —  + 

m 


1 


1  +  £(1  -coajd 


(4) 


k2n(k.  H) 

(i  +  ^(i-cos^))2 

We  study  Llie  case  when  ^  ^  1.  i,e,  the  photon  energy  is  much  less  then 
the  electron  rest  energy.  In  the  first  approximation  over  the  ratio  of  these  two 
energies  the  photon  distribution  function  does  not  change  in  the  collision  as 
the  function  of  the  wave  number  value  k*  It  changes  only  as  tile  function  of 
the  scattering  angle.  Thus  (4)  could  l>e  integrated  over  k .  Finally  the  kinetic 
equation  have  the  form: 


dn  \/l  —  p2  d(rn)  dn 
dt  ^  r  dr  +  ^  dz 


=  'RtlN(z)  {(3 /.yJ  -  n2  - ^  +  3)(n(/i')  -  nM))pl  •  (5) 

where  ji  —  C0$i 9.  //  =  CQ$$\  iS-  -  polar  anges  of  the  income  and  outcome 
photons  respectively. 

Function  n(t,  r.  Zyfj)  is  expanded  over  Legendre  polinoms  in  order  to  simplify 
the  collision  integral  of  photons  with  electrons 

n{L.r.z,p)  =  y~Vg(t,r,  c)Pfe(/r) 
fe=0 
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After  substitution  of  these  expressions  into  (5)  and  elementary  transformations 
the  sought  kinetic  equation  could  Ire  obtained: 


dn  yfi  —  /i2  d(rn)  dn 
dt  r  dr  ^  dz 


(fi) 


=  (TThN(z) 


no  -ft  +  —  n2P2(/i) 


)• 


I  fere  JV(^)  -  the  air  density,  07^  the  cross-section  of  the  Tbmsou  scattering  of 
the  photons  on  electrons,  /i  -  the  cosine  of  the  angle  between  the  photon  impulse 
direction  .and  the  iion-uiiifomity  gradient,  -  photon  distribution  function 
averaged  over  the  angles,  n  2-  photon  distribution  function  inter  gated  together 
with  the  second  Legendre  polinoni,  -  second  Legendre  polinom.  Note, 

that  no  ass  options  about  tlie  value  of  the  non- uniformity  gradient  was  done 
while  obtaning  tlie  kinetic  equation.  So,  the  equation  (6)  as  well  as  the  initial 
equation  (1)  is  true  for  the  arbitrary  value  of  atmosphere  noii-unifomity.  The 
analysis  of  the  equation  ff>)  was  started.  It  was  shown  that  its  exact  solution  for 
tlie  function  t)  is  completely  determined  only  by  the  zero  and  second 

expansion  coefficients  of  the  function  over  Legendre  polinoms. 

After  substitution  of  these  expressions  into  (5)  and  elementary  transformations  the  sought 
kinetic  equation  should  be  obtained: 


dn  <yjl-  jU2  d{rn) 


dn 


~  +  V^  =  °ThNtZ)\ 

or  oz 


nn  —  n  H - 

0  10 


(6) 


dt  r 

Here  the  air  density  changes  with  the  altitude  as 

N(z)  =  N0e  a\ 

where  No  -  the  electron  density  at  the  height  z  —  0,  Va  -  the  characteristic  scale  of 


atmosphere  density  change.  n  -  the  cross-section  of  the  Tomson  scattering  of  the  photons  on 
electrons,  tomcohobckoc  ceLieimc  paccomioi  (Jiotohob  xieicipoiiaMii.  //  -  the  cosine  of  the 

angle  between  the  photon  impulse  and  the  non-uniformity  gradient,  ^  (/0  -  second  Legendre 
polynom.  Note,  that  no  assumptions  about  the  value  of  the  non-uniformity  gradient  was  done 
while  obtaining  the  kinetic  equation.  So,  the  equation  (6)  as  well  as  the  initial  equation  (1)  is 
true  for  the  arbitrary  value  of  atmosphere  non-uniformity.  The  exact  solution  for  the  function 

n(z,  r,  //,  t)  js  completely  determined  only  by  the  zero  and  second  expansion  coefficients  of 
the  function  over  Legendre  polynoms. 

Introducing  a  new  mark  for  °ThN o,  equation  (6)  could  be  transferred  to  the  form: 


dn  dri  1  az  1 

—  +  ii—  +  7e  n  =  -re  az 
ot  oz  l  l 

Consider  the  linear  differential  operator 


(7lO  +  £712^2  (/*))• 


(7) 


Ln{t,  z,n)=  (  —  +  n—  +  -c~aZ  J  AO  =  N(t’  z>  tl)‘ 

The  solution  of  this  inhomogeneous  partial  differential  equation  could  be  obtained  by  the 
standard  Green  function  method. 

t  H-co  1 

n(t,z,fi)  =  n(0,  z,n)  4-  J  j  J  d^dzfd^fG(tfz,  ,zr,f/), 

o  -oo-l 
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where  n(0,  z,  /J)  -  n(t,  z,  fi)  at  the  initial  moment  (t,  —  0),  G{t-  z,  jr.  t' ,  z1 ,  //) .  Green  function 
obeying  the  following  equation 

{§t  +  fijL  +  le~asS)  G  =  s&-  “  *')*(**  ~ 

(  eaz  \ 

Presenting  function  G  as  a  product  exP  )9  one  can  obtain 

(l + *4) 9 = exp  (~^r)s(t  ■ m 3  ■ 

Now  transforming  9  as  exP  where  ^  obey  the  equation 

§t  + 4  =  5{t  -  e)5{*  -z’^-  <9  > 

Its  solution  has  the  form 

X  =  Q(t  —  t')5(z  —  zf  —  /i(t  —  tf))d(^  —  {/). 

It  means  that  the  Green  function  is  equal  to 

C(t,  z ,  //;  y,  /)  =  exp  (c~°^°~OZ  )  X 


x5(z  —  z'  —  /r(<  — 1/))5(//  —  —  0> 

Thus,  the  differential  equation  (7)  could  be  rewritten  as  the  integral  equation 

t,  H-no  -T  T  *  _  _  f  -v 

n(t, z,n)  =  n(0T z,/u)  +  7/  /  f  dt'dz'dy! exp  fB  “  -azrJx 

0  -oo  -1  V  7 

X$(z  -  z'  -  fl(t  -  t'))5(fl  -  fl')  [no{tf,z’)  +E7l2(t',  Z,)P2(flf)]  . 

The  analysis  of  equation  (9)  allowed  to  find  the  form  of  the  photon  distribution  function  on 

condition  that  ^ !mc  >>  1 .  where  k-  photon  momentum.  The  distribution  function  is  found  as 
a  function  of  an  angle,  time  and  height.  The  form  the  function  for  different  values  of  parameter 
az  and  time  moments  at  is  shown  in  Fig.  3.2. 


(8) 


(9) 


Fig.3.2.  /(A  zi  I-1)  as  the  function  of  !l  at:  a)  az  —  10,  at  =  1,  al  —  0.1 ;  b)  az  —  10,  at  =  1, 

al  —  0.01;  c)  az  —  5,  at  =  1,  al  —  0.1. 

The  angle  dependence  of  the  photon  density  n(t,z,fi )  is  determined  by  the  function 

6(ti)  =  p  +  n2, 
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(2_M)~28 


where 

This  implies  that  the  preferential  direction  in  the  photon  propagation  iz  along  the  axis  of 
inhomogenity  Oz  towards  the  less  medium  density  (Fig.3.3). 


t* 

Fig.3.3.  Plot  of  W. 
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Task  4:  High-energy  particles  in  thunder  (lightning)  phenomena.  Experimental  and 
theoretical  investigation  of  intensive  gamma  emission,  generation  of  HF  and  VHF  radio 
waves,  electric  field,  electric  field  variations 


Experiment 

DISCRETION  OF  THE  THUNDERSTORM  COMPLEX 
The  experimental  complex  "Thunderstorm"  is  located  at  the  height  3340  m  above  the  sea 
level.  It  consists  of  the  following  installations:  storm  trigger  system,  scintillation  system  of 
Nal-detectors,  multi-layers  spectrometer  of  absorption,  muon  hodoscopes,  monitoring  system 
of  high-energy  and  thermal  neutrons,  two  independent  radio  systems  working  in  a  range  of 
frequencies  0,1  -  30  MHz  and  ~  250  MHz,  detector  of  jump  of  a  static  electrical  field  and  its 
high-frequency  component.  All  experimental  complex  works  continuously  from  the  moment  of 
start  (April)  before  ending  a  season  of  thunderstorms  (September)  in  an  automatic  mode. 

The  storm  trigger  system  consists  of  several  tens  detectors  created  on  the  basis  of  Geiger 
counters  SI5G.  This  system  fixes  the  moment  of  extensive  atmosphere  showers  of  cosmic  rays 
(EAS)  passage,  allows  estimating  its  size,  and  also  energy  of  a  primary  particle.  One  EAS 
detector  consists  of  20  counters  SI5G,  and  its  sensitive  area  is  about  0,6  m  .  The  EAS  trigger 
detectors  are  covered  with  iron  filter  10  cm  width  (Fig. 4. 2)  in  order  to  exclude  accidental 
coincidences  induced  by  low-energy  (below  10  MeV)  electrons. 


Fig.4.1.  Relative  location  of  EAS  triggers  registration  points. 
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Fig. 4.2.  EAS  trigger  system  detectors. 

The  EAS  trigger  is  produced  when  the  pulses  from  four  neighbor  registration  points 
coincide  in  time.  The  scheme  of  the  trigger  formation  is  presented  in  Fig.4.3. 

The  computer  modeling  showed  that  such  configuration  of  EAS  trigger  system  provides 
the  effective  selection  of  EAS  having  primary  energy  higher  than  1  PeV  and  which  axes  cross 
the  system  at  the  distance  up  to  100m. 

The  scintillation  system  from  14  scintillation  Nal-detectors  is  used  for  registration  of 
intensity  of  soft  y  and  x-ray  radiation  from  thunderstorm  clouds.  Detectors  are  situated  in  pairs 
in  7  registration  points  build  on  the  Tyan-Shan  Station  pass  and  on  the  flanks  of  surrounding 
hills.  They  are  placed  as  a  chain  (see  Fig.4.3)  perpendicular  to  the  usual  direction  of  the 
thundercloud  movement.  The  distance  between  the  extreme  points  in  this  chain  is  about  2  km, 
the  vertical  interval  occupied  by  the  registration  points  is  about  500  m.  Thus  the 
“Thunderstorm”  scintillation  system  allows  studying  directly  the  space  distribution  of  the 
emission  inside  thunderclouds  both  in  horizontal  and  vertical  directions. 
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Fig.4.3.  Location  of  the  Nal  scintillation  detectors  on  the  Tyan-Shan  Station  pass. 


Two  scintillation  detectors  are  situated  in  each  registration  points.  Nal  crystal  in  one  of 
this  detector  is  placed  inside  the  metal  (aluminum,  1  mm)  case.  The  case  of  the  second  detector 
is  done  of  a  bit  of  the  polyethylene  pipe  with  wall  thickness  10mm.  As  a  result  two  detectors 
have  different  gamma  emission  registrations  thresholds.  In  addition  in  each  registration  points 
most  distant  from  the  center  (PUP1,  PUP2,  TUR2  h  TUR3  on  the  Figure)  one  “empty”, 
without  Nal  crystal,  detector  is  placed.  The  signals  of  “empty”  detectors  are  linked  to  the  joint 
registration  system.  They  serve  for  the  control  of  electromagnetic  breakthroughs  on  cables 
connecting  the  registration  points  with  the  data  control  center. 

Pulses  with  standard  amplitude  (TTL  compatible)  formed  by  comparators  are  5  ps  long.  They 
are  transferred  to  the  data  registration  system  based  on  CAM  AC.  The  system  provides 
recording  of  signal  intensity  during  4000  time  intervals  having  duration  200  ps  each.  This 
system  operates  continuously.  It  keeps  in  memory  the  intensity  time  dependence  for  last  0.8  s. 
When  the  special  control  signal,  the  trigger,  come  the  data  registration  system  continues  its 
work  for  0.4  s.  After  that  the  data  accumulated  are  recorded  on  the  disk.  Thus,  the  time 
dependence  of  signal  intensity  is  recorded  with  200  ps  resolution  during  0.4  s  after  the  trigger 
moment. 

The  absorption  spectrometer  consists  of  layers  SI5G  counters  alternated  by  thin  layers 
Pb  and  Fe.  Absorption  spectrometer  is  used  as  for  the  registration  of  electrons  accelerated  in 
thunderstorm  electric  field  as  for  the  registration  of  y-emission  from  these  electrons.  The 
energy  of  electrons  and  y-quanta  is  estimated  as  well. 

«THUNDERSTORM»  complex  includes  «Radio-HF»  installations  for  radio  emission 
waveform  measurements  in  the  frequency  range  of  0,1  -  30  MHz  and  «Radio-E»  installation  to  measure 
electric  field  and  its  variations. 

«Radio-HF»  installation  is  designed  to  record  short  (down  to  30  ns)  electromagnetic  pulses 
generated  in  thunderclouds  by  lightning  discharges  and  their  sources  direction  finding.  It  consists  of 
three  spaced  (by  about  70  m)  antenna  assemblies,  connected  to  central  recording  unit  by  cables  of  equal 
length.  Each  antenna  assembly  contains  three  antennas:  two  crossed  by  90°  loops  to  measure  horizontal 
magnetic  field  component,  and  End-Fed  antenna  to  measure  vertical  component  of  electric  field. 

The  complex  includes  two  such  installations  spaced  by  about  400  m  that  provides  a  possibility  to 
determine  radio  emission  source  location  at  distances  up  to  2  -  4  km.  Radio  installations  operate  in  the 
external  trigger  mode.  The  overall  record  duration  was  200  ms  with  prehistory  time  160  ms  (record  time 
duration  before  the  trigger  coming). 
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«Radio-E»  installation  is  designed  to  record  electric  field  under  thunderstorm  conditions,  its 
variations  is  the  frequency  range  0.5-25  kHz  and  250  MHz  radio  emissions.  It  also  produces  the  drive 
pulse  when  the  amplitude  of  electric  field  variation  overrides  the  threshold  value.  Static  (“slow”) 
electric  field  is  measured  by  electrostatic  flux  meter  ("field  mill"  type).  Electric  field  variations  («fast» 
field)  are  measured  using  capacitor  type  sensor.  At  the  lightning  discharge  moment  «Radio-E» 
equipment  registers  static  electric  field  change,  its  fast  variations  (stepped  leader  and  return  stroke  for 
cloud-to-ground  discharges),  and  form  a  control  signal  (lightning  trigger)  for  other  installations  of  the 
experimental  complex.  The  installation  uninterruptedly  registers  data  in  the  “slow”  mode  with  the 
frequency  of  20  Hz.  After  the  first  drive  pulse  is  produced  it  turns  into  the  “fast”  registration  mode.  In 
this  mode  the  data  digitization  rate  is  20  jus  during  one  hour,  after  that  the  installation  turns  back  to  the 
“slow”  mode. 

Data  processing 

The  database  containing  events  registered  at  the  Tien-Shan  Mountain  Station  during  2007 
summer  observation  campaign  was  created.  It  is  based  on  the  MySQL  5.1  SQL-server  and  use 
OpenOffice.org  2.4  as  a  SQL-client.  For  each  registered  event  the  following  data  were  put  in: 

1.  date; 

2.  time; 

3.  registration  of  short  wave  radio  emission  (SW); 

4.  registration  of  electric  field  jump; 

5.  registration  of  gamma  emission; 

6.  registration  of  gamma  emission  spectrum; 

7.  the  type  of  the  trigger  (EAS,  fast  electric  field  jump,  electromagnetic  pulse)  initiated 
the  record; 

8.  the  estimation  of  the  electric  field  jump; 

9.  qualitative  estimation  of  the  gamma  emission  (weak  or  strong); 

10.  registration  of  the  noise  in  control  gamma  registration  channels; 

11.  registration  of  SW  emission  by  each  of  Radio- 1  and  Radio-2  installations; 

12.  registration  of  the  electric  field  jump; 

13.  registration  of  gamma  emission  in  each  of  the  observation  points  TUR1,  TUR2,  TUR3, 
KAP,  PUP1,  PUP2,  NEV. 

All  the  events  registered  during  the  thunderstorm  activity  on  08/01,  08/02,  08/03,  08/04,  08/07, 
08/08  and  08/25,  2007  are  placed  in  the  database. 


Registered  events,  total 

169 

Including: 

0 

Caused  by  the  electromagnetic  pulse  trigger 

968 

Caused  by  the  EAS  trigger 

686 

Caused  by  the  fast  electric  field  jump  trigger 

36 

The  small  number  of  fast  electric  field  jump  trigger  events  is  determined  by  the  high  threshold 
for  this  trigger  fixed  in  2007. 

The  selection  of  the  similar  events  for  the  purpose  of  the  further  analysis  was  done  using 
database.  Numbers  of  events  in  the  selected  groups  are  presented  in  the  Table. 


Electromagnetic  pulse  trigger  +  SW  emission  is  registered 

601 

EAS  trigger  +  gamma  emission  is  registered 

623 

EAS  trigger  +  gamma  emission  is  registered  +  spectrum  of  gamma  emission  is 

49 
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registered 

EAS  trigger  +  SW  emission  is  registered 

36 

EAS  trigger  +  gamma  emission  is  registered  +  SW  emission  is  registered 

29 

Gamma  emission  is  registered  +  electric  field  jump  is  registered 

192 

Gamma  emission  is  registered  +  gamma  emission  is  registered  +  SW  emission  is 
registered 

131 

Strong  gamma  emission  is  registered  +  SW  emission  is  registered 

115 

Gamma  emission  is  registered  in  two  or  more  registration  points 

685 

Strong  gamma  emission  is  registered  in  two  or  more  registration  points 

512 

The  Section  THUNDER  of  the  site  http://www.tien-shan.org/she/vardbaccess/  was 
developed  (in  Russian).  This  Section  contains  the  time  scans  of  radio  emission  registered  in 
2007  in  all  registration  points  as  well  as  the  observation  log.  Fast  view  of  the  gamma  emission 
at  registration  points  TUR1,  TUR2  h  TUR3  is  available. 

The  data  processing  showed  up  the  short  gamma  flashes  events  having  complex 
structure:  10-30  MeV  electrons;  low  energy  (dozens  keV)  and  high  energy  (~MeV)  gamma 
quanta.  There  is  a  temporary  correlation  of  short-term  flashes  with  the  moments  of  the 
lightning  registered  by  electromagnetic  pulse  trigger.  The  spatial  correlation  of  these  flashes  of 
radiation  with  the  presence  of  thunderclouds  nearby  the  experimental  setup  location  is  found. 
The  spatial  sizes  of  area  of  flare  are  about  hundred  meters.  Probably  this  type  of  events  is 
connected  with  cloud-to-cloud  electric  breakdown. 


35 


i  i  j _  T7-|naj  rT^»r*V^nirja1  (t?  ppnr(t 


^6 


P 

:: 

1 1 

£ 

CC. 

itt 

=> 


£ 


p  ri  50  ir'Kju  1950  iuOG  2050  2iOO  2iSO 

NUMin.n  OF-  t  fMiu  i Nj  1 1 i-ivAL 


£ 

q: 

m 
I  L 


y* 


fMI  06  2007  12:06  57  Mu6  at  TUR1 


IB50  1900  1950  2000  2050  2100  2150 


NUMBER  OF  TIME  INTERVA1 


Fig.4.4.  Gamma  emission  bursts  registered  by  scintillate  counters  triggered  by  lightning. 
Intensities  are  shown  for  three  registration  points  (TUR-III  (540  m  above  the  Tien-Shan  Station 
level);  TUR-II  (310  m);  KAPT  (Tien-Shan  Station  level).  The  intensities  for  each  counter  are 
shown  for  three  registration  thresholds,  which  correspond  to  the  minimum  energy  of  registered 
quanta  of  15,  60  and  240  keV. 
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Electron  fluxes 

During  the  2009  thunder  season  on  Tien-Shan  Cosmic  Ray  Station  the  emission  intensity 
was  measured  using  the  multi-layer  absorption  spectrometers  bases  on  SI5G  ionization 
counters.  Unlike  the  measurements  done  in  1999-2004  the  present  ones  were  done  using  the 
system  of  time  base  which  allowed  microsecond  resolution  registering  of  the  intensity  change 
with  time.  The  three-layer  spectrometers  in  detector  points  NEV  and  100M  were  used. 

The  ionization  counters  SI5G  were  placed  in  each  of  the  detector  points  in  three  layers 
(“up”,  “middle”  and  “lower”)  situated  one  above  the  other.  Each  layer  consists  of  60  counters. 
The  sensitive  area  of  each  counter  (n  x  D  x  L)  is  0.1 1  m2,  the  total  sensitive  area  of  one  layer  is 
about  6  m2.  Counters  are  situated  in  the  wood  boxes  done  out  of  the  5-mm  thick  plywood.  A 
30  mm  thick  rubber  absorber  separates  the  upper  and  middle  layers  while  the  middle  and  lower 
layers  are  separated  with  1  mm  lead  sheet.  The  same  sheet  covers  the  floor  of  the  detector  point 
screening  from  below  the  whole  detector  assemblage  from  the  base  emission.  The  output 
signals  of  all  counters  of  the  same  layer  are  bridged  and  are  transferred  to  the  registering 
system  through  the  common  cable. 

Thus  3  signals  are  received  by  the  data  acquisition  system  from  each  of  detector  points 
NEV  n  100M.  They  are:  pulses  from  upper,  middle  and  lower  layers  of  ionization  counters. 

Two  additional  signals,  the  coincidence  (within  5  ps  time  gate)  pulse  between  the  counter 
layers  (“upper”  +  ’’middle”  and  ’’upper”  +  ’’middle”  +  “lower”)  are  formed  in  the  data 
acquisition  system  itself.  The  efficiency  of  gamma  quanta  registration  in  SI5G  counters  is  0.05 

-  2%  depending  on  the  energy  while  their  efficiency  of  charged  particle  registering  is  about  95 

-  99%.  That’s  way  the  registering  of  the  coincidence  signals  allows  to  distinguish  the  pulse 
from  charged  particles  (electrons)  from  those  from  gamma  quanta  in  the  total  pulse  flow. 

The  data  acquisition  system  allows  obtaining  the  time  scans  of  the  intensity  during  4000 
serial  time  intervals  having  duration  200  ps  each.  The  moment  of  the  record  of  the  scan  is 
determined  by  a  special  signal  -  trigger  which  position  in  the  time  scan  corresponds  to  the  time 
interval  having  number  2000.  Thus  the  time  dependence  of  the  intensity  is  represented  during 
2000  time  intervals  (0.4  s)  before  trigger  moment  and  2000  time  intervals  after  this  moment. 

The  trigger  was  formed  in  the  same  way  as  the  electromagnetic  trigger  in  2007  year 
measurements.  The  trigger  was  formed  as  coincidence  pulse  of  electromagnetic  pulses  induced 
on  four  long  (500-1000  m)  cables  extended  along  the  Station.  At  the  end  of  the  season  (during 
thunderstorms  on  /08/30/09  and  /09/01/09)  the  shower  trigger  system  was  used  as  well.  It 
formed  the  trigger  when  the  extensive  atmosphere  showers  crossed  simultaneously  its  six 
widespread  detector  points. 

The  first  type  of  registered  events  (see  for  example  Fig.4.5)  is  presented  by  time  scans  in 
wich  short  (having  duration  1  -10  time  intervals)  growths  of  intensity  in  lower  layer  (some 
time  even  in  middle)  is  seen  while  there  was  no  any  noticeable  signal  in  the  upper  layer  and  no 
any  coincidence  signals.  These  events  should  be  interpreted  as  the  registration  of  the 
bremsstrahlung.  The  electrons  accelerated  in  the  thunderstorm  electric  field  remain  not 
registered  in  the  upper  cause  the  Compton  effect  and  the  photo  effect  in  the  absorber  situated 
lower.  The  electrons  arising  in  the  process  are  registered  by  the  counters  situated  in  lower 
layers.  At  that  the  energy  of  gamma  emission  converts  more  effectively  in  the  lead  sheet 
determining  the  intensive  signal  from  counters  in  the  lower  layer. 

The  cross  section  of  photo  effect  rapidly  decreases  with  the  energy  of  electromagnetic 
emission.  That’s  way  the  existence  of  such  kind  of  events  means  that  gamma  quanta  with 
energies  up  to  1  MeV  are  registered  in  the  experiment.  It  corresponds  to  characteristic  energies 
of  bremsstrahlung  from  accelerated  electrons  and  confirms  our  results  of  2007  received  by  the 
use  of  Nal  detectors. 

Such  type  of  events  form  the  main  part  of  registered  events,  they  are  observed  exclusively 
in  the  active  stage  of  a  thunderstorm  while  they  are  completely  missed  during  fine  weather  as 
well  as  during  the  rains  and  snowfalls  (even  intensive)  without  thunderstorms. 
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The  signals  from  electromagnetic  emission  are  divided  into  two  groups  according  to  their 
durations:  short  flashes  having  duration  1-5  time  intervals  (i.e.  less  then  1  ms)  and  relatively 
long  flashes  having  duration  more  then  10  time  intervals,  i.e.  several  microseconds.  This 
distinction  corresponds  to  that  observed  in  2007  in  scans  of  gamma  intensity  measured  by  the 
use  of  Nal  scintillators:  short  intensive  flashes  of  gamma  emission  originated  from  the 
bremsstrahlung  of  electron  avalanches  and  more  long  flashes  of  softer  gamma  emission 
registered  usually  in  the  events  caused  by  the  electromagnetic  trigger. 

Electromagnetic  flashes  overlap  the  moment  of  the  trigger  in  about  half  part  of  all  events. 
They  come  upon  the  time  intervals  with  numbers  ±10  -  20.  It  means  that  these  flashes 
correspond  to  the  moments  of  the  maximum  electric  field  strength  which  occur  immediately 
before  the  lightning  generating  the  trigger. 

In  the  second  type  of  events  (see  Fig. 4. 6)  the  narrow  spikes  of  intensity  (with  characteristic 
duration  -  one  time  interval,  i.e.  200  ps)  are  observed.  These  spikes  simultaneously  are 
displayed  in  the  signals  from  upper,  middle  and  lower  layers  of  ionization  counters  as  well  as 
in  the  coincidence  signals  of  both  types.  The  probability  to  register  series  of  quanta  of 
electromagnetic  radiation  in  two  or  even  in  three  counter  layers  is  negligible.  Thus,  the 
existence  of  coincidence  signals  means  that  these  spikes  correspond  to  cases  of  registration  of 
electrons  accelerated  in  the  field  of  a  thundercloud. 

The  events  of  registering  of  charged  particles  occur  much  often  then  the  events  with  spikes 
in  time  scans.  As  opposite  to  the  electromagnetic  spikes  scans  with  charged  particles  are 
observed  only  in  one  detector  point.  Thus,  the  characteristic  path  of  electrons  between  these 
poins  is  about  150- 160m.  The  events  with  charged  particles  just  as  events  of  the  first  type  are 
observed  exclusively  in  the  active  stage  of  a  thunderstorm  while  they  are  completely  missed 
during  fine  weather  as  well  as  during  the  rains  and  snowfalls  (even  intensive)  without 
thunderstorms. 

The  base  level  of  intensity  for  the  coincidence  signals  is  0.07  -  0.1  pulses  per  one  interval 
of  the  time  scan.  The  typical  number  of  coincidence  signals  is  3-5  pulses  per  interval,  it  for  30- 
70  times  exceeds  the  base  level.  The  placement  of  these  signals  on  the  time  scale  is  rather 
uniform  what  confirms  the  year  2007  observations  when  narrow  intensity  spikes  did  not 
concentrate  to  the  trigger  moment  also  (in  contradictory  to  the  long  flashes  of  soft  emission 
with  the  events  with  electromagnetic  trigger). 

During  the  thunderstorm  on  September  1,  2009  both  electromagnetic  trigger  and  shower 
trigger  were  used.  The  last  formed  the  trigger  when  the  extensive  atmosphere  shower  (EAS) 
passed  through  six  detectors  (A,  B,  C,  D,  E,  F)  widespread  along  the  Station.  The  trigger  was 
formed  in  the  case  of  coincidence  of  at  least  four  signals  from  detector  points  A-F  within  2  ps 
time  gate. 

In  some  time  scans  of  shower  trigger  events  obtained  during  the  thunderstorm  on 
September  1,  2009  the  intensive  flashes  of  electromagnetic  emission  and  the  flows  of  charged 
particle  correlated  with  these  flashes.  The  flows  of  charged  particle  are  observed  nearby  the 
moments  of  trigger  left  behind  it  for  1-2  time  intervals.  Most  likely  these  intensive  flashes  are 
observed  in  the  cases  when  the  EAS  direction  was  favorable  and  the  flow  of  particles 
registered  by  the  shower  trigger  system  passed  through  the  thundercloud.  Events  of  such  kind 
could  be  considered  as  a  direct  observation  of  the  lightning  initiation  by  the  EAS  particles. 
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Fig.4.7 

Theory 


The  theory  of  the  lightning  stepped  leader  which  is  under  development  in  the  Project 
conjointly  takes  into  account  three  main  processes:  conventional  breakdown,  phenomenon  of 
electron  runaway  and  runaway  breakdown. 

STEPPED  LEADER 

The  lightning  leader  steps  have  widely  varied  characteristics.  According  to  the  statements  in 
the  textbooks  [1-2]  based  on  the  numerous  measurements  one  can  indicate  the  following  main 
features  of  the  steps. 

(1)  The  step  length  is  3-50  m.  As  example  one  can  take  an  average  value  15  m.  Note  that  the 
length  is  slightly  growing  with  the  atmospheric  height. 

(2)  The  step  pulse  duration  is  of  the  order  0.3-1  ps. 

(3)  The  time  between  pulses  either  or  between  steps  changes  randomly  from  5  to  100  ps. 

(4)  The  pulse  velocity  obtained  from  the  measurements  of  the  luminosity  movement  is  ~5x  109 
cm/s. 

(5)  The  width  of  the  luminous  region  is  up  to  1  m. 

(6)  The  radio  emission  measured  at  the  Earth  surface  demonstrate  that  the  leader  step  pulse  is 
accompanied  by  a  significant  pulse  of  electric  current  /  ~  0.3-5  kA. 


CONVENTIONAL  BREAKDOWN 

Let  us  determine  the  basic  pulse  characteristics  from  the  model  of  electric  discharge  in  a  strong 
electric  field  E  ~  \0Eth  ~  200  kV/cm.  The  conventional  breakdown  in  a  given  strong  electric 
field  leads  to  a  very  fast  exponential  growth  of  air  ionization  Ne 

Ne  =  N0  exp (viont)  (1) 
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Here  No  is  the  initial  density  of  electrons  and  vi0„  is  the  ionization  frequency.  The  breakdown  in 
air  was  studied  in  details  in  numerous  experimental  and  theoretical  works  (see  textbooks  [4-5]). 
The  ionization  frequency  is  growing  effectively  with  the  electric  field  value: 

Vion  =  VamE(E  /  E»h  )'  (2) 

Here  vam  is  the  maximal  attachment  frequency  in  e-C>2  collisions: 


v,_  =7.6x10“13AL 


-i 


2x10 


yFmO  J 


-1 


Nm0  =2.7x10 


19  cm  3 


A  well-known  empirical  formula  [4]  for  F  is 

■ 5/- 

F  =  I  E/ 


(3) 


(4) 


It  is  correct  at  low  values  of  E/Eth  up  to  E/Eth  ~4-5  (see  [16]).  At  higher  values  of  (E/Eth)  the 
growth  of  vion/vam  is  diminishing  reaching  maximum  ~  1 05  approximately  at  (E/Eth)  ~  100.  In 
the  region  we  are  interested  in 


=7-10, 


=  (2-5)xl03 


(5) 


So,  we  see  from  (l)-(5),  that  due  to  the  breakdown  in  a  given  strong  electric  field  the 
degree  of  ionization  in  air  grows  very  rapidly. 

On  the  other  hand,  the  growth  of  the  free  electron  density  leads  to  the  fast  amplification 
of  electron  conductivity  o.  The  last  process  leads,  as  it  is  well  known,  to  the  diminishing  of 
electric  field: 

E  =  E0  exp(-  4^|  odt^  (g) 

which  will  stop  the  growth  of  the  ionization.  Electrical  conductivity  in  neutral  atmosphere 


a  = 


e2Ne 

m  v 


(7) 


where  vem  is  the  electron  collision  frequency.  Substituting  (6),  (7)  in  (1)  we  find  that  electric 
field  keeps  it’s  high  value  for  the  finite  time  t  <  t*  when  N(t)  <  N(t*): 

47re2N(t*) 

(8) 
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This  relation  determines  the  maximal  value  of  electron  density  reaching  in  breakdown: 


N 


F(- 
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4  ne“  /  ^  th 

and  characteristic  time  of  the  discharge 
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AT 


N. 


■  10-20. 


(9) 


(10) 


Relations  (9),  (10)  determine  the  basic  characteristics  of  a  strong  electric  field  (E»Eth) 
discharge  in  unbounded  air 


RUNAWAY  ELECTRONS  IN  HIGH  ELECTRIC  FIELD 

The  electron  runaway  effect  in  gas  under  the  action  of  the  electric  field  is  connected  with  the 
decreasing  of  the  breaking  force  with  the  increasing  of  the  electron  energy  s'. 

F  a  s<jrNZ  oc  — 

K  a 


(11) 
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Here  Na  -  atomic  number  density,  Z  -  charge  of  a  nuclear,  ^  R 


<7 R  >  \hr  Rutherford  cross- 
section.  In  the  strong  electric  field  E  an  electron  of  sufficient  high  energy 
eE 


£>£0  = 


(JnZN, 


(12) 


will  be  continuously  accelerated.  Such  electrons  are  called  runaway  electrons.  If  the  electric 
field  is  not  too  strong  E«Ecn  the  main  part  of  electrons  having  low  energy  e  <  £,  (where  £,-  is 
the  ionization  energy)  are  colliding  with  neutral  atom  as  a  whole.  Friction  force  is 


(13) 


F=°aNaVT’ 

where  m  and  va  electron  mass  and  collision  frequency 

Va=(7aNaVT ,  (14) 

vT  =  yfTjm — electron  thermal  velocity,  oa — collision  cross  section.  Eq.  (13)  determines  usual 
drift  of  electrons  under  the  action  of  electric  field  E  =  F/e.  Thus  we  see  that  at  low  electron 
energies  s  <  the  friction  force  of  electrons  moving  in  the  neutral  gas  F  is  growing  with  their 
energy.  But  at  high  energies  s  >  sm  it  begins  to  fall  down.  It  means  that  at  some  middle 
energies  s  ~  sm  ~Zsi  it  will  reach  the  maximal  value.  We  emphasize  that  around  this  region  the 
interaction  of  fast  electron  with  neutrals  gradually  transforms  from  elastic  and  inelastic 
scattering  at  the  atom  as  a  whole  to  the  Coulomb  scattering  at  the  electrons  and  the  atomic 
nucleus.  To  describe  this  process  it  is  necessary  to  use  the  kinetic  theory. 

Kinetic  theory  of  runaway  electrons  in  neutral  gasses  was  developed  by  Gurevich  [6]. 
The  electric  field  E  was  supposed  homogeneous  and  the  gas  unbounded.  According  to  the 
theory  the  maximal  value  of  friction  force  Fmax  is  characterized  by  the  electric  field  Ecn: 


4m3NZ 


(15) 


(16) 


cn  e  2.72/ 

where  1  is  an  average  excitation  energy.  Taking  into  account  that  I~  13. 5Z  eV,  an 
astonishingly  simple  expression  for  the  value  of  critical  field  in  neutral  gas  was  obtained: 

Ecn  =7  xlO~15 Na  V/cm. 

Here  Na  is  the  number  density  of  atoms  in  cm3.  For  example  in  air 

iV  ) 

Ecn  =38[27xi0i9  j  MV/m  *20 Eth 
(17) 

where  Nm  is  the  number  density  of  air  molecules,  and  E,h — a  threshold  electric  field  of 
conventional  breakdown. 

If  the  electric  field  is  less  than  critical  one  E  <  Ecn  only  the  electrons  in  the  tail  of 
distribution  function  become  runaways.  The  flux  of  runaway  electrons  is  proportional  to  the 
exponential  factor  [6]: 

<i8> 

Here  Ne  is  the  number  density  of  electrons,  ve — electron  collision  frequency  at  a  characteristic 
runaway  electron  energy  ec  and  A  is  a  large  constant  (for  air  A  ~  30).  It  means  that  the  number 
of  runaway  electrons  falls  down  very  rapidly  with  electric  field  decreasing.  The  basic  theory 
was  farther  developed  and  confirmed  for  different  gases  in  a  number  theoretical  and 
experimental  studies  (see  textbooks  [7,8]). 


RUNAWAY  BREAKDOWN  IN  STRONG  ELECTRIC  FIELD 


The  basic  physical  process  determining  the  phenomenon  of  runaway  breakdown  (RB)  is  the 
generation  of  new  fast  electrons  due  to  the  runaway  particle  ionization  of  neutral  molecules. 
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Although  the  majority  of  newborn  free  electrons  have  low  energies,  some  will  have  rather  high 
energy  s  >ec.  As  a  result,  an  exponentially  growing  runaway  avalanche  can  occur.  Two  main 
features  of  RB  are: 

(1)  It  can  occur  in  the  electric  field  E  >EC,  where  Ec  is  an  order  of  magnitude  less  then  the 
threshold  field  of  conventional  breakdown  Ec  «Eth. 

(2)  The  presence  of  “seed”  fast  electrons  having  energies  above  the  critical  runaway  energy  sc 
is  needed  to  initiate  RB. 

In  the  previous  works  the  RB  theory  in  weak  electric  fields  E  «  Eth  was  constructed  and  a 
role  of  RB  in  lightning  leader  initiation  was  studied  [10-11],  But  one  can  show  that  the 
runaway  breakdown  effect  can  take  place  in  a  strong  electric  field  E  »  Eth  as  well.  Our  goal  is 
to  develop  such  a  theory  -  the  theory  of  strong  runaway  breakdown  (SRB).  It  will  be  shown 
that  usual  runaway  electrons  serve  as  a  seed  for  SRB,  when  their  energy  reach  the  runaway 
critical  energy  ec  which  coincide  exactly  with  the  one  determined  in  the  weak  RB  theory: 
Sc=(mc2/2\EJE ). 

Thus  the  usual  runaway  effect  in  neutral  gases  gradually  transforms  to 
SRB,  at  electron  energies  s: 

£>£c=  (me2  2\EC  E) «  25^*  /  j  keV  (19) 


Here  we  took  into  account  that  for  air  Ec/Eth  ~  0.1.  It  is  evident  from  (17)  and  (18)  that  SRB 
effect  at  E  ~  (l-4)Eth  is  extremely  low.  It  explains  that  in  usual  long  sparks  E  <  (3-4)  Eth  (see 
[12])  the  runaway  electrons  and  gamma  emission  are  practically  absent.  On  the  contrary,  at 
E/Eth>l-\0  the  runaways  could  become  significant  what  lead  to  the  exponentially  strong  SRB. 
Temporal  and  spatial  growth  of  runaway  breakdown  in  strong  electric  field 
Let  us  consider  the  strong  runaway  breakdown  (SRB)  effect  in  a  strong  electric  field  E/Et h  ~  7- 
10.  It  means  that  parameter  d  =  E/Ec  ~  70-100.  In  such  strong  fields  we  can  consider  the  main 
part  of  distribution  function  in  nonrelativistic  limit.  Moreover,  the  distribution  function  is 
directional  along  electric  field.  The  kinetic  equation  describing  runaway  breakdown  in  this  case 
has  a  form  [13] 


dt 


1  df 

-^r—  + 

v2  dv 


1 

Av5 


J  f(k)kdk. 


(20) 


Here  z — direction  along  electric  field,  A — Bethe’s  logarithm.  Eq.  (20)  has  a  very  important 
feature:  the  parameter  §  could  be  eliminated  from  equation.  Actually,  introducing 

u 

F  =  \f(k)kdk,  (21) 


X  -  zS2 ,  T  =  tS2 

we  can  rewrite  Eq.  (20)  in  the  form 

d2F  1  d2F  Id 
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u  =  vS 


dvdx  u  dudr  u2  du 


yU  du  y 
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This  self-similarity  of  Eq.  (23)  on  parameter  S  means  that  we  can  solve  it  for  d  =  1  and 
automatically  find  from  (22)  the  runaway  breakdown  growth  rates  for  arbitrary  S. 

Let  us  consider  first  the  limit  case  of  spatial  uniform  time  growth. 

In  this  case  we  can  search  for  the  solution  of  Eq.  (20)  in  the  form / (v,  t)  =f  (v)  exp(7/rrb4 
According  to  (21),  (23)  the  time  of  exponential  growth  of  electrons 

Trb  =  T0<r3/2. 

This  result  is  in  agreement  with  spatially  uniform  runaway  breakdown  theory  [13]. 

(2)  Spatial  dependence 

In  this  case  we  can  search  for  the  solution  of  Eq.  (20)  in  the  form / (v,  t)  =f  (v)  exp(z//a). 
According  to  (21),  (23)  the  avalanche  scale  4  decrease  with  increasing  S 

4  =  IW2 


(22) 

(23) 


(24) 


(25). 


44 


3641p  Final  Technical  Report 

It  means  that  the  avalanche  length  in  a  strong  field  is  diminishing  dramatically.  Thus,  in  air  for 
S  =  70-100  we  have  4  ~  1  cm.  The  distribution  function  for  S  =  70  and  the  growth  of  runaway 
electron  density  for  S  =  70,  100  determined  by  numerical  solution  of  Eq.  (23)  are  presented  at 
Figs.  4.8,  4.9.  The  results  are  in  agreement  with  studied  previously  runaway  breakdown  for  the 
case  of  a  weak  electric  field.  The  distribution  function  is  approaching  rapidly  to  f~  1/e.  The 
characteristic  time-scale  of  SRB  process  in  atmosphere  is  of  the  order  of  a  few  ns.  Note  that  the 
density  growth  is  stronger  during  the  initial  0.5  ns,  analogous  effect  was  seen  in  the  weak  fields 
(see  [13]). 


Fig.4.8.  Dependence  of  distribution  function  of  fast  electrons  on  energy  for  S  =  100.  Different 
curves  correspond  to  time  moments  0.01,  0.02,  0.03  and  0.04.  Asymptotic  behavior  f~  1  /e  is 
clearly  seen. 
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Fig.4.9.  Growth  of  total  number  of  fast  electrons  with  time  for  parameters  8  =  100  and  8  =  70. 
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Numerical  simulations  of  runaway  breakdown  in  inhomogeneous  medium 

The  electron  generation  during  runaway  breakdown  in  the  case  of  spatial  non-uniformity  was 
studied.  The  computations  are  performed  by  means  of  the  program  for  solving  the  kinetic 
equation  that  takes  into  account  the  dependence  of  the  electron  distribution  function  on 
coordinate.  It  is  found  that  the  electron  density  increases  sharply  near  the  point  where  the 
electric  field  is  maximal.  The  region  of  the  electron  density  increasing  is  small  enough.  Here 
the  flow  velocity  of  electrons  is  abruptly  increasing.  The  spatial  distributions  of  the  electron 
density  and  flow  velocity  are  illustrated  in  Fig.4.10.  The  distributions  are  given  at  time  t= 0.05 
(dashed  and  dot-dashed  curves,  respectively)  and  t=0.1  (solid  and  dotted  curves,  respectively). 
The  electric  field  amplitude  E/Eo=5  0. 

Numerical  computations  show  that  the  initial  electron  flow  velocity  has  a  small  effect 
on  the  rate  of  the  electron  generation.  The  effect  depends  on  the  electric  field  amplitude.  Also, 
there  is  a  small  effect  of  the  size  of  the  electric  field  region. 


Fig.4.10.  Spatial  distributions  of  the  electron  density  and  flow  velocity 

The  numerical  investigation  of  the  runaway  breakdown  on  the  front  of  the  relativistic  wave  was 
studied.  The  increment  of  the  particle  number  increasing  strongly  increase  with  the  increasing 
of  the  electric  field  amplitude  E/Eo,  where  iso  -  the  critical  field  of  breakdown.  It  is  clearly  seen 
in  Fig  4.1 1  were  the  time  dependences  of  the  electron  density  are  presented  for  different  values 
of  E/E0.  Since  the  time  t=0.01  electron  density  exponentially  increases  with  the  constant 
increment.  Its  value  is  in  a  good  agreement  with  the  theoretically  expected  dependence  0.35 
(E/Eof  2.  The  deviations  to  the  more  small  values  are  observed  at  relatively  small  values  of 
E/Eo<l5.  Electrons  with  small  energies  y< 3  give  the  main  contribution  into  the  density. 
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Fig.4.1 1.  Growth  of  fast  electron  number  as  a  function  of  the  electric  field  amplitude. 


Task  5:  Transformation  of  an  extensive  atmospheric  shower  in  the  thunder  electric  field: 
intensive  growth  of  electromagnetic  component,  namely  of  the  high-energy  electrons 
number,  of  gamma  quanta  and  positrons  number.  Space  enlargement  of  the  shower. 
Experimental  and  theoretical  investigation 


Experiment 

During  a  summer  season  thunderstorm  clouds  pass  directly  at  the  height  of  the  Station 
(3340  above  see  level),  so  the  installation  detectors  are  located  deeply  inside  of  a  cloud.  In  the 
summer  season  2007Y  the  on-line  monitoring  of  extensive  air  showers  was  carried  out  and 
their  correlations  with  flares  of  gamma  and  X-radiation  were  studied.  The  goal  of  the 
investigations  was  to  check  theoretical  predictions  of  combined  effects  of  runaway  breakdown 
and  high-energy  cosmic  rays  ( E  >  1015-1016  eV)  on  lightning  processes  in  the  thunderstorm 
atmosphere. 

The  measurements  with  scintillation  Nal-detectors  included  registration  of  the  intensity 
in  slow  and  fast  (100-200  ps)  time  span.  Measurements  with  the  high  resolution  were 
conducted  in  4000  consecutive  time  slots  having  200  ps  duration;  2000  intervals  before  and 
2000  after  the  trigger  moment.  The  measurements  were  carried  out  using  the  system  of  Nal- 
detectors  with  the  wide  site  separation  (-2  km  in  a  horizontal  and  -500-600  m  in  a  vertical 
direction).  Short  flares  of  radiations  were  found  for  the  first  time  in  measurements  with  the 
high  temporary  resolutions.  They  have  complex  structure  indicating  the  presence  of  electrons 
accelerated  till  10-30  MeV  as  well  as  electrons  having  low  (tens  keV)  and  high  (deciles  and 
units  of  MeV)  energies. 


3641p 


48 


Final  Technical  Report 

Three  types  of  trigger  signals  were  used  in  measurements  with  the  high  resolution:  the 
EAS  trigger,  ’’jump  of  a  field”  trigger  and  electromagnetic  pulse  trigger.  Four  separate  systems 
of  concurrences  constructed  on  the  base  hodoscopes  from  SI5G  Geiger  counters  were  used  for 
a  trigger  pulse  creation  at  the  moment  of  EAS  passage  through  the  installation.  The  EAS 
trigger  pulse  initiated  the  data  recording  from  all  subsystems  of  a  ’’THUNDERSTORM” 
complex. 

The  lightning  discharge  electromagnetic  emission  was  studied  using  two  radio  installations  in 
the  0.1-30  MHz  frequency  range.  They  allowed  registering  the  pulse  waveform  with  16  ns 
time  accuracy  as  well  as  to  determine  the  direction  of  the  emission  source.  The  installation 
“Radio-E”  was  used  in  the  measurements  as  well.  The  operation  threshold  was  chosen  too 
large  in  the  2007  season,  thus  it  corresponded  to  the  nearby  discharges  only.  That  why  the  total 
number  of  the  produced  drive  pulses  was  small  compared  to  the  previous  seasons  when  the 
threshold  was  noticeably  lower  and  false  events  where  registered  due  to  the  noise  even  without 
thunderstorm.  Radio  installations  “Radio-HF”  were  set  going  by  the  trigger  pulses  from  the 
“Radio-E”  installation  as  well  as  from  the  storm  trigger  system.  Moreover,  installations  could 
be  set  going  spontaneously  due  to  the  interference  induction  from  the  lightning  discharges, 
most  likely  on  the  flux,  -  electromagnetic  pulses  (EMP).  It  could  be  supposed  that  at  least  in 
some  cases  these  are  the  same  pulses  that  set  going  the  scintillate  detector  record  system.  It 
must  be  mentioned  that  not  all  of  the  trigger  pulses  could  be  processed  by  the  radio 
installations.  It  concerns  with  the  large  data  amount  obtained  in  the  single  event.  Data  must  be 
recorded  of  the  hard  disk  of  the  computer,  the  record  of  one  event  takes  30  seconds,  at  that  the 
record  time  varies  from  one  vent  to  the  other  and  varies  for  different  installations.  The  trigger 
pulses  coming  during  the  recording  are  not  processed.  The  comparison  of  data  obtained  in 
simultaneous  measurements  of  gamma  and  radio  emission  showed  that  as  a  rule  both  emissions 
are  observed  in  those  cases  when  it  could  be  supposed  that  the  startup  of  corresponding 
installations  takes  place  due  to  the  same  EMP.  A  record  example  with  the  assumption  of  the 
synchronous  operating  time  of  installations  is  shown  in  Fig.5.1. 
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The  high  time  correlation  between  gamma  and  radio  emission  is  seen  in  the  event.  In 
Fig. 5.2  quasi-stationary  electric  field  in  the  same  time  range  is  shown.  The  minor  field  jump  is 
observed  corresponding  to  the  distant  cloud-to  ground  discharge.  It  should  be  mentioned  that 
the  time  lock-on  of  the  data  obtained  on  the  “Radio-E”  installation  to  the  other  installations 
could  be  done  in  2007  campaign  with  not  best  than  the  second  accuracy. 
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Fig.  5.2. 

Gamma  and  radio  emission  is  observed  in  events  corresponding  to  the  EAS  trigger 
rarely  than  in  the  events  corresponding  to  the  EMP  trigger,  radio  emission  is  observed  more 
often  than  gamma  one.  It  is  determined  by  the  fact  that  the  passage  of  the  shower  is  not 
specially  or  temporally  connected  with  the  existence  of  strong  electric  field  in  the 
thundercloud,  typical  for  the  lightning  discharge  initiation.  According  to  that  most  showers  on 
the  one  hand  can’t  initiate  the  runaway  breakdown  effect  in  a  cloud  and  on  the  other  hand  can’t 
be  affected  by  the  thunder  electric  field.  In  the  opposite  case,  when  the  shower  passes  through 
the  strong  electric  field  region  of  the  thundercloud  which  exceeds  the  threshold  electric  field  of 
runaway  breakdown,  the  EAS  could  be  the  initiator  of  the  lightning  discharge  generating 
powerful  E-field  and  gamma  emissions.  It  could  be  observed  at  that  if  sensors  are  not  far  from 
the  region  where  EAS  interacts  with  the  cloud,  due  to  the  small  gamma  ray  path  length  in 
atmosphere.  In  opposite,  the  radio  emission  could  be  observed  from  the  large  distance  from  the 
interaction  region.  An  example  of  the  probable  lightning  initiation  by  the  EAS  is  shown  in  Fig. 
5.3,  the  discharge  beginning  coincides  in  time  with  the  shower  passage.  The  gamma  emission 
in  this  event  was  not  observed,  most  likely  due  to  its  remoteness  from  the  sensor  placement. 
An  example  of  event  corresponding  to  the  EAS  trigger  is  shown  in  Fig. 5. 4.  Both  radio  and 
gamma  emission  are  observed,  they  are  well  correlated.  It  should  be  mentioned  that  in  this 
event  the  registered  shower  did  not  initiate  the  lightning  discharge,  the  emission  started  earlier, 
but  it  have  through  the  region  of  the  strong  electric  field  where  the  discharge  was  developed.  It 
appeared  in  the  existence  of  intensive  pulse  of  gamma  and  radio  emission  at  the  moment 
corresponding  to  the  shower  passage  time. 
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The  system  of  gamma-radiation  detectors  permits  to  analyze  the  amplitude  information 
of  the  input  scintillation  signals.  To  use  this  information  the  absolute  detector  calibration  is 
needed.  A  connection  between  the  six  thresholds  of  pulse  discriminators  of  the  registration 
system,  and  corresponding  ranges  of  energy  deposit  of  gamma-radiation  quanta  inside  the 
scintillating  crystal  was  determined.  After  calibration  it  was  established  that  pulse 
discriminators  have  monotonously  increasing  operation  thresholds,  corresponding  to  gamma- 
quanta  energy  range:  30,  40,  60,  70,  120  and  320  keV. 

The  spatial  distribution  of  bursts  observed  in  our  experiment  is  about  1  km  as  it  is  seen 
simultaneously  at  four  detector  points  for  all  energy  thresholds.  We  have  compared  the  data  for 
different  energy  thresholds  at  different  detector  points  and  found  that  the  integral  energy 
spectrum  is  similar.  It  allowed  us  to  construct  the  averaged  integral  spectrum  taking  into 
account  the  results  of  detector  calibration.  The  spectrum  is  presented  at  the  Fig.5.5.  It 
demonstrates  a  good  agreement  of  our  experiments  both  with  the  Eack  et  al.  balloon 
experiments  [5]  and  with  runaway  breakdown  theory. 


Figure  5.5.  Integral  energy  spectrum  of  short-time  gamma  flashes  in  the  thunderstorm 
atmosphere:  data  of  Tien-Shan  experimental  (squares),  theoretical  curve  (solid  line)  and 
balloon  experiment  data  by  Eack  et  al.  [5]  (triangles). 

The  EAS  triggered  gamma  bursts  were  studied.  The  typical  example  of  such  an  event  is 
presented  at  Fig  5.1.  The  burst  of  gamma  emission  of  the  type  as  in  Fig  5.6  time  scans  lasts  for 
about  one  millisecond.  This  burst  is  the  powerful  single  outburst  of  hard  gamma  quanta 
registered  in  all  three  energy  diapasons.  Besides,  time  scans  in  Fig  5.6  demonstrate  the  wide 
space  distribution  of  the  emission  of  the  burst.  The  signal  is  observed  in  registration  points 
TUR3,  TUR2  and  KAPT  without  significant  depression.  It  disappears  only  in  lower  point 
NEV. 

The  probability  to  observe  the  EAS  triggered  burst  in  the  time  scan  increases  from  2  % 
in  fair  weather  to  10%  when  thundercloud  do  exist.  The  total  fraction  of  bursts  registered  in 
time  scans  is  quit  a  few  (only  10%  even  in  the  thunder  case)  among  all  registered  EAS.  It  could 
be  apparently  explained  by  the  relatively  small  area  of  the  shower  trigger  system  used  in  the 
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experiment  and  its  large  distance  from  those  registering  points  where  the  most  powerful  signals 
from  bursts  were  observed,  i.e.  from  “upper”  scintillators  TUR3,  TUR2  and  TUR1. 

The  temporary  correlation  of  short-term  flares  with  the  moments  of  EAS  passage 
through  area  with  high  intensity  of  an  electrical  field  inside  clouds  is  found.  The  flashes 
observable  at  the  EAS  trigger  have  much  smaller  duration  (parts  of  ms)  compared  to  the  case 
of  the  electromagnetic  trigger  and  are  formed  by  gammas  with  ~  MeV  energy.  Spatial  area  of 
these  flashes  is  much  more  wide  ~  1  -2  km. 
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Fig. 5. 6.  EAS  triggered  event:  time  scans  of  the  gamma  emission  intensity  at  the  height 
540  m  (TUR3)  and  310  m  ({TUR2)  above  the  average  level  of  the  Tien-Shan  Station.  Abscissa 
axis:  the  number  of  the  time  interval,  200  ps  each. 
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08.03.2007  12:25:16  Mu4  at  KAPT+EAS 
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Phc.5.6.  (Continued).  EAS  triggered  event:  time  scans  of  the  gamma  emission  intensity  at  the 
height  60  m  (KAPT)  and  at  the  Tien-Shan  Station  level  (NEV). 
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During  11  hours  of  thunderstorm  in  August,  2007  we  observed  600  events  triggered  by  EAS. 
Strongly  amplified  pulse  of  gamma  radiation  simultaneously  with  EAS  was  seen  in  4  events. 
One  event  happened  15.08.07  during  initial  phase  of  thunderstorm,  other  three  events  took 
place  08.08.07  during  the  main  phase  of  the  storm.  Two  events  will  describe  here:  one  in  the 
initial,  another  one  in  the  main  phase. 

Event  15.08.07 

Storm  lasted  15.08.07  for  four  hours  from  05.18  up  to  09.40.  The  event  took  place  at  05:18:12 
in  the  initial  phase  of  the  storm.  The  shower  was  fixed  by  all  six  detectors  of  EAS  trigger 
system.  That  means  the  EAS  was  generated  by  cosmic  ray  particle  with  the  energy  not  less 
than  1016  eV.  The  pulse  of  gamma  emission  was  seen  at  the  scans  of  three  observational  points: 
7,  6  and  4.  The  scans  are  presented  in  the  Fig. 5. 7. 


15.08.2007  -  05:18:12 
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Fig.5.7.  Time  scans  of  gamma  emission  for  the  event  15.08.07.  Each  time  interval  is  equal  to 
200  mcs.  The  moment  of  the  EAS  trigger  corresponds  to  2000-th  time  interval. 
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08.08.2007-  12:13:35 
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Fig.  5.8.  Time  scans  for  the  event  08.08.07.  From  up  to  down:  gamma  emission  measured  by 
detectors  ScI-5,  ScI-6  and  ScII-7.  The  bottom  scan:  F1F  radio  emission. 

The  zero  is  the  moment  of  the  EAS  trigger. 
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One  can  see  from  the  Fig.  5.7: 

•  The  strong  pulse  of  gamma  emission  was  observed  in  all  points  simultaneously 

•  The  pulse  width  is  one  bin,  less  than  200  mcs 

•  The  moment  of  gamma  pulse  in  all  observational  points  coincide  exactly  with  EAS 

•  A  strong  change  of  gamma  ray  background  was  observed. 

Note,  that  the  distance  between  observation  detector  points  in  space  is  about  1 1 50m. 
They  are  distributed  in  height  at  600m.  So,  the  strong  pulse  of  gamma  emission  covered  quite  a 
large  area.  One  can  see  also  from  the  Fig.5.7  that  the  observed  energy  distribution  of  gamma 
quanta  has  a  classical  RB  form:  small  change  between  40  keV  (1st  thr.)  and  60  KeV  (2st  thr) 
and  strong  diminishing  to  120  KeV  (3st  thr). 

Of  special  interest  is  the  sharp  change  of  the  gamma  background,  seen  at  the  scan  in 
Fig.5.7.  One  can  see  that  the  background  value  at  the  point  6  for  three  times  diminishes  at  the 
trigger  moment  exactly.  At  the  point  7  the  background  is  growing  up  during  50  mcs  before  the 
trigger  moment.  Then  it  falls  down  approximately  for  three  times.  The  growth  of  the  gamma 
breakdown  produced  by  cosmic  rays  is  a  consequence  of  the  increasing  of  the  electron  flux  in 
the  thunderstorm  electric  field  due  to  RB  process.  The  sharp  fall  of  the  background  means  that 
the  electric  field  of  the  thundercloud  crossed  by  EAS  has  sharply  diminished.  In  other  words, 
an  electric  discharge  of  the  cloud  took  place.  A  combined  action  of  EAS  and  RB  gave  a 
macroscopic  effect.  So,  this  is  just  the  observation  of  RB-EAS  discharge. 

Event  8.08.07 

The  main  phase  of  a  storm  8.08.07  with  frequent  lightning  discharges  lasted  from  1 1:53 
to  12:52.  The  event  under  discussion  was  registered  at  12:13:35.  The  shower  was  fixed  by  four 
detectors  of  the  EAS  trigger  system  what  means  that  the  energy  of  cosmic  ray  particle  was 
about  1015  eV.  The  pulse  of  gamma  radiation  was  measured  by  detectors  at  the  observational 
points  5,  6  and  7. 

The  scans  are  shown  in  the  Fig.5.8.  The  strong  gamma  pulse  simultaneous  with  EAS 
trigger  at  all  observational  points  is  seen.  The  measured  radio  signal  is  also  shown  in  the 
Fig.5.8.  One  can  see  a  very  strong  pulse  of  radio  emission  at  the  EAS  trigger  moment  exactly. 
It  means  that  the  interaction  of  EAS  with  the  cloud  was  accompanied  by  the  electric  discharge. 

Discussion 

The  probability  of  RB-EAS  event  observation  is  determined  by  two  factors: 

•  The  trajectory  of  cosmic  ray  particle  which  generate  RB-EAS  event  should  cross  the 
thundercloud  at  the  distance  400-500  m  from  gamma  detectors.  Otherwise  gamma 
mission  cannot  be  observed  due  to  the  absorption  in  the  atmosphere.  On  the  other  hand 
the  same  trajectory  should  come  to  our  EAS  trigger  system  to  be  fixed  as  a  shower. 
These  conditions  limit  the  quantity  of  the  possible  trajectories. 

•  The  second  condition  is  natural:  the  value  of  electric  field  in  a  thundercloud  crossed  by 
EAS  should  be  higher  than  the  RB  critical  field. 

This  condition  is  evidently  more  often  fulfilled  during  the  main  phase  of  the  storm. 

Due  to  these  two  factors  RB-EAS  is  a  rare  event.  It  is  of  the  order  of  1%  from  all  EAS 
events,  observed  during  the  thunderstorm.  This  estimate  agrees  with  the  results  of  our 
measurements:  4  RB-EAS  events  among  600  EAS.  Thus  the  simultaneous  observation  of  EAS 
with  a  strong  gamma  emission  pulse  in  a  wide  space  region  and  probability  considerations 
demonstrates  the  experimental  verification  of  the  RB-EAS  phenomenon,  predicted  for  a 
thunderstorm  atmosphere. 

Note,  that  in  a  fine  weather  conditions  analogous  strong  bursts  of  gamma  emission  with 
EAS  were  never  registered. 
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The  presented  features  of  short-term  flares  allow  to  state  that  they  arise  during  the 
development  of  an  electron  cascade  in  a  thundercloud  electric  field  with  the  subsequent 
slowdown  of  accelerated  electrons  in  an  atmosphere,  i.e.  they  are  the  direct  experimental 
confirmation  of  the  realization  of  the  RB-  effect. 


Theory 

Note  on  theory  of  runaway  breakdown  (RB)  effect 

Asymptotic  solution 

The  problem  of  matter  breakdown  is  formulated  in  the  kinetic  theory  in  the  following 
form.  The  kinetic  equation  for  electrons  is  considered  in  the  coordinate  and  momentum  space 
(r,  p).  The  electric  field  E,  the  neutral  molecular  density  N,  their  charge  Z  and  other  substance 
parameters  determining  the  electron  collisions  are  considered  as  constants  in  the  simplest 


problem  statement.  Collision  integral  is  linearized.  Then  asymptotic  over  time  t  solution  for  the 
distribution  function  /( PT)  of  the  kinetic  equation  has  the  exponential  form: 

/( P,  0  /oi  <P)  exp(v )  +  /02  (p)  exp(v2t) 


It  is  uniform  and  does  not  depend  on  space  coordinates.  Here  vi  -  eigenvalue  of  the  linear 

kinetic  equation,  i= 1,2,  foi  (P)  represents  the  corresponding  eigensolution. 

The  exponential  growth  of  the  distribution  function  and  consequently  the  growth  of  the 

electron  number  just  mean  the  breakdown  of  the  substance.  Parameter  vi  determines  the 

ionization  frequency  and  vi  ~  Ti  -  the  characteristic  time  of  breakdown.  The  existence  of  two 
independent  solutions  of  the  linear  kinetic  equation  means  that  there  are  two  types  of 
breakdown  in  any  dielectric:  the  conventional  breakdown  (R)  and  RB  (R). 


Similarity  relations 

The  runaway  breakdown  can  take  place  in  any  substance.  The  interaction  of  fast 
electrons  is  coulomb  and  consequently  it  always  has  a  similar  character.  Thus  RB  has  the  same 
structure  for  all  substances  and  possesses  the  remarkable  similarity  properties.  Namely,  the 
critical  field  Ec  is  proportional  to  the  mass  density  p  in  all  substances.  If  p  is  expressed  in 
grams  per  cm3,  then 

Ec  =l,8/?[MeFcm  '] 

Correspondingly,  the  critical  energy  of  runaway  s  is  connected  with  the  value  of  the 
electric  field  E  by  the  relation 

me2 

s  =  — - — , 

28 

where  ^  =  E/ Ec  _ 

The  characteristic  length  of  breakdown  l  strongly  decreases  with  the  growth  of  the 
electric  field  (~  $  2) ,  while  the  ionization  frequency  correspondingly  increases: 

T2=\Q-wp-l8~2l2[c\ 

T,=W°p52l2[c-x) 

In  dense  substances  with  /?~10-100  g  cm-3  the  characteristic  times  of  the  breakdown  are 
extremely  small.  It  should  be  added  that  these  times  rapidly  fall  down  with  the  growth  of  the 
parameter  <£ 

Note,  that  the  indicated  similarity  relations  are  correct  in  the  limited  range  of  variation 
of  the  parameters): 

1,5  <  <5  <  100-150. 
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Runaway  breakdown  in  atmosphere 

The  dependence  of  RB  ionization  frequency  17  for  air  on  the  electric  field  value  E  is 
presented  in  Fig.5.3.  One  can  see  that  RB  ionization  frequency  monotonically  increases  with 
the  field  E  (~  Eyi ).  At  the  same  time  the  conventional  breakdown  ionization  frequency  17  at 
first  increases  very  rapidly  with  the  increasing  of  E  (approximately  as  ~  E5’5 ),  but  then 
saturates  (see  [1-3]). 

It  is  important  that  RB  can  take  place  even  at  low  values  of  electric  field,  i.e.  in  the  small  field, 
when  the  condition  Eth  >  E  >  Ec  js  fulfilled.  In  the  case  the  ionization  frequency  is  extremely 
high:  v2  >  107  c’1. 

The  conventional  breakdown  can  not  take  place  in  the  small  field,  it  occurs  at  E  >  Eth . 
Moreover,  although  the  ionization  frequency  of  the  conventional  breakdown  very  rapidly 
increases  with  the  increasing  E  it  still  remains  less  than  the  RB  ionization  frequency  up  to  the 

value  E  -  2Eth  This  circumstance  can  be  highly  important  and  can  lead  to  the  occurrence  of 

the  high  number  of  fast  electrons  in  the  conventional  discharge  at  ^ Elh  >E>Eth  p  may  be 
assumed  that  it  is  the  cause  of  the  observed  bursts  of  gamma  ray  emission  both  in  a  lightning 
leader  [4-6]  and  in  discharges  studied  in  laboratories  [7-9], 

So,  RB  is  the  main  breakdown  in  air  relative  to  the  quickness  of  ionization  increasing 

not  only  at  E  <  Eth  but  at  larger  values  of  electric  field  up  to  E  =  2 Eth 

The  second  important  parameter  presented  in  Fig.5.9  is  the  characteristic  length  /.  It 
determines  the  minimal  space  scale  of  the  breakdown  range.  One  can  see  that  in  small  fields 
the  length  is  very  large  (/  «  30  -  50  m).  But  it  rapidly  decreases  with  the  increasing  of  electric 
field  E  (~  E~2),  so  that  even  for  E  ~  (1  -  2 )Eth  the  length  /  becomes  of  the  order  of  10-30  cm. 

At  high  field  values  E  »  Eth  the  characteristic  length  /  becomes  rather  small.  It  demonstrates 
that  while  at  small  electric  fields  the  investigation  of  RB  is  possible  only  in  thunderclouds,  the 
investigation  at  higher  electric  fields  is  possible  on  laboratory  installations  as  well. 
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Fig.5.9.  Dependence  of  the  conventional  breakdown  ionization  frequency  17  and  of  the  RB 
ionization  frequency  V2  on  the  electric  field  value  (left  scale,  c_1). 
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Dependence  of  the  avalanche  growth  length  l  on  the  electric  field  value  (right  scale,  m)  is 

presented  as  well. 

Effect  of  cosmic  rays 

To  RB  occur  besides  the  fulfillment  of  the  condition  E  >  Ec  the  presence  of  fast  seed 

electrons  having  energy  exceeding  the  runaway  energy  ( 8  >  e c )  is  needed.  These  conditions 
are  realized  in  thunderclouds.  Observations  show  that  maximal  electric  fields  in  thunderclouds 
are  close  to  the  critical  RB  field  Ec  [1,2,9],  The  secondary  electrons  of  cosmic  rays  (CR)  serve 

as  fast  seed  particles.  Their  mean  flux  density  at  heights  4-8  km  is  relatively  high:  (E<  ~  103 . 

The  essential  difference  of  RB  from  conventional  breakdown  consists  in  the  effective 
generation  of  X  and  gamma  emission.  Exactly  the  observation  of  gamma  emission  with 
s  —  50  1 00  keV  indicates  the  possible  RB  realization.  Simultaneously  strong  lectric  currents 

could  be  excited,  leading  to  the  generation  of  radio  emission. 

It  is  natural  that  the  number  of  runaway  electrons  is  proportional  to  the  number  of  seed 
ones  generated  in  the  interaction  of  CR  with  the  air.  The  total  seed  electron  number  Ne  in  an 
EAS  increases  proportional  to  the  energy  of  the  primary  CR  particle.  For  example,  Ne  is  about 

106  in  the  shower  generated  by  the  primary  CR  particle  with  the  energy  ECR  «  1018  eV.  If 

£c**1018ev,then^»1010. 

While  an  EAS  passes  through  the  thundercloud  where  E  >  Ec  the  runaway  electrons 
avalanche  grows  up  and  the  number  of  energetic  electrons  in  the  EAS  increases  exponentially. 
Simultaneously  the  number  of  thermal  electrons  increases  for  several  millions  times.  All 
together  they  generate  the  RB-EAS  discharge  [10].  Naturally,  this  discharge  should  be 
accompanied  by  the  strong  pulse  of  gamma  emission  due  to  the  large  number  of  fast  electrons. 
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